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We show that, as a result of non-linear self-interactions, it is feasible, at least in light of the bounds 
coming from terrestrial tests of gravity, measurements of the Casimir force and those constraints 
' imposed by the physics of compact objects, big-bang nucleosynthesis and measurements of the 

| cosmic microwave background, for there to exist, in our Universe, one or more scalar fields that 

couple to matter much more strongly than gravity does. These scalar fields behave like chameleons: 
£SJ ■ changing their properties to fit their surroundings. As a result these scalar fields can be not only very 

strongly coupled to matter, but also remain relatively light over solar system scales. These fields 
could also be detected by a number of future experiments provided they are properly designed to 

Qdo so. These results open up an altogether new window, which might lead to a completely different 
view of the role played by light scalar fields in particle physics and cosmology. 
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I. INTRODUCTION 



There is wide-spread interest in the possibility that, in addition to the matter described by the standard model 
of particle physics, our Universe may be populated by one or more scalar fields. These are a general feature in 
high energy physics beyond the standard model and are often related to the presence of extra-dimensions, [lj. The 
existence of scalar fields has also been postulated as means to explain the early and late time acceleration of the 
Universe 0, 0, 0, B IS 0] • It is almost always the case that such fields interact with matter: either due to a direct 
Lagrangian coupling or indirectly through a coupling to the Ricci scalar or as the result of quantum loop corrections, 
[1, H, [lO, EH, E2|- If the scalar field self-interactions are negligible, then the experimental bounds on such a field 
are ver y st rong: requiring it to either couple to matter much more weakly than gravity does, or to be very heavy 
[II EI El Ell- Recently, a novel scenario was presented by Khoury and Weltman [17( that employed self-interactions 
of the scalar-field to avoid the most restrictive of the current bounds. In the models that they proposed, a scalar field 
couples to matter with gravitational strength, in harmony with general expectations from string theory whilst, at the 
same time, remaining very light on cosmological scales. In this paper we will go much further and show, contrary to 
most expectations, that the scenario presented in [l7| allows scalar fields, which are very light on cosmological scales, 
to couple to matter much more strongly than gravity does, and yet still satisfy all of the current experimental and 
observational constraints. 

The cosmological value of such a field evolves over Hubble time-scales and could potentially cause the late-time 
acceleration of our Universe [Hj]. The crucial feature that these models possess is that the mass of the scalar field 
depends on the local background matter density. On Earth, where the density is some 10 30 times higher than the 
cosmological background, the Compton wavelength of the field is sufficiently small as to satisfy all existing tests of 
gravity. In the solar system, where the density is several orders of magnitude smaller, the Compton wavelength of 
the field can be much larger. This means that, in those models, it is possible for the scalar field to have a mass in 
the solar system that is much smaller than was previously thought allowed. In the cosmos, the field is lighter still 
and its energy density evolves slowly over cosmological time-scales and it could function as an effective cosmological 
constant. While the idea of a density-dependent mass term is not new H3, HH, HH, [23|, [24j, the work presented 
in [It], EH is novel in that the scalar field can couple directly to matter with gravitational strength. If a scalar field 
theory contains a mechanism by which the scalar field can obtain a mass that is greater in high-density regions than 
in sparse ones, we deem it to possess a chameleon mechanism and be a chameleon field theory. When referring to 
chameleon theories, it is common to refer to the scalar field as the chameleon. 

We start this article by reviewing the main features of scalar field theories with a chameleon mechanism. Afterwards, 
this paper is divided into roughly two parts: in sections IIII1 IIVI and fVl we study the behaviour of chameleon theories as 
field theories, and derive some important results. From sections I VII onwards, we combine these results with a number 
of experimental and astrophysical limits to constrain the unknown parameters of these chameleon theories {(3, M, A}. 
We shall show how the non-linear effects, identified in sections IllUfVl allow for a very large matter coupling, (3, to be 
compatible with all the available data. We also note that some laboratory based tests of gravity need to be redesigned, 
if they are to be able to detect the chameleon. If the design of these experiments can be adjusted in the required way, 
and their current precision maintained of its current level, then a large range of sub-Planckian chameleon theories 
could be detected, or ruled out, in the near future. 

In section [TTTl we study how <fi behaves both inside and outside an isolated body and derive the conditions that 
must hold for such a body to have a thin-shell. In this section, we show how non-linear effects ensure that the value 
that the chameleon takes far away from a body with a thin-shell is independent of the matter-coupling, [3. Whilst 
such (3- independence as been noted before for </> 4 -theory in [37[, this is the first time that it has been shown to be a 
generic prediction of a large class of chameleon theories. In section HVl we show the internal, microscopic, structure of 
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macroscopic bodies can unexpectedly alter the macroscopic behaviour of the chameleon. Using the results of section 
IIIII and IIVI we are then able to calculate the </>-force between two bodies; this is done in section [V] In each of these 
sections we take care to note precisely when linearisation of the chameleon field equation is invalid. 

Laboratory bounds on chameleon field theories are analysed in section l*VTl We focus mainly on the Eot-Wash exper- 
iment reported in [38l 139 ] , which tests for corrections to the 1 /r 2 behaviour of gravity, and experimental programmes 
that measurement the Casimir force [13, EH, |42| • We also look at the variety of laboratory and solar-system based tests 
for violations of the weak equivalence principle (WEP) [26|, [62|, [63|, [64[ . The extent to which proposed satellite-based 
searches for WEP violation will aide in the search for scalar fields with a chameleon-like behaviour is considered in 
this section. We shall see that for a large range of values of M and A, laboratory tests of gravity at unable to place 
any upper-bound on (i. 

In sections IVIII and I Villi we show how the stability of white-dwarfs and neutron stars, as well as requirements 
coming from big bang nucleosynthesis and the Cosmic Microwave Background, can be used to bound the parameters 
of chameleon field theories. We shall see that such considerations do result in upper-bounds on j3. 

Finally, in sections Hxl and Ixl we collate all of the different experimental and astrophysical restrictions on chameleon 
theories, use them to plot the allowed values of /3, M and A, and discuss our results and their implications. 

We include a summary of the main results at the end of sections IlllfVl for easy reference. This allows the reader, 
less interested on the detailed derivation of the formulae, to follow the whole article. Throughout this work we take 
the signature of spacetime to be (H ) and set H = c = 1; G = 1/M?j where M p i is the Planck mass. 

II. CHAMELEON FIELD THEORIES 

In the theories proposed in [T3|, the chameleon mechanism was realised by giving the scalar field both a potential, 
V(<f)), and a coupling to matter, B{(5(f)/M p i)p, where p is the local density of matter. We shall say more about how 
the functions V and B are defined, and the meaning of /3, below. The potential and the coupling-to-matter combine 
to create an effective potential for the chameleon field: V e ^ {(f) — V((f>) + B{(3cf)/M p i)p. The values <fi takes at the 
minima of this effective potential will generally depend on the local density of matter. If at a minima of V e ^ we have 
(j) = 4> c , i.e. V e ^ (4> c ) = 0, then the effective 'mass' (m c ) of small perturbations about <f> c , will be given by the second 

derivative of V e ^ ', i.e. to 2 = V^((f> c ). It is usually the case that |V^| ^> \B^p\ and so m c will be determined 
almost entirely by the form of V{4>) and the value of <p c . If V(<p) is neither constant, linear nor quadratic in (f> then 
V^(<p c ), and hence the mass m c , will depend on (f> c . Since 4> c depends on the background density of matter, the 
effective mass will also be density-dependent. Such a form for V(cj>) inevitably results in non-linear field equations for 

For a scalar field theory to be a chameleon theory, the effective mass of the scalar must increase as the background 
density increases. This implies V^tf, (4>c) /Vfiific) > 0. It is important to note that it is not necessary for either 
V(cf>), or B((3(p/Mpi), to have any minima themselves for the effective potential, V e '\ to have minimum. A sketch of 
the chameleon mechanism, as described above, is shown in FIGS. [T]and[2j In FIG [1] the potential is taken to be of 
runaway form and has no minimum itself. However, It is clear from the sketches that V e ^ does have an minimum, 
and that the value <f> takes at that minimum is density dependent. In FIG [5] the potential is taken to behave like 4 
and so does have a minimum at (f> = 0. However, the minimum of the effective potential, V e ^ , does not coincide 
with that of V . Once again, the minimum of V e ^ is seen to be density dependent. 

A. The Thin-Shell Effect 

This chameleon mechanism often results in macroscopic bodies developing what is called a "thin-shell" . A body is 
said to have a thin-shell if </> is approximately constant everywhere inside the body apart from in a small region near 
the surface of the body. Large (0(1)) changes in the value of 4> can and do occur in this surface layer or thin- shell. 
Inside a body with a thin-shell V</> vanishes everywhere apart from in a thin surface layer. Since the force mediated 
by <j) is proportional to V</>, it is only that surface layer, or thin-shell, that both feels and contributes to the 'fifth 
force' mediated by <j>. 

It was noted in [T3, EH that the existence of such a thin-shell effect allows scalar field theories with a chameleon 
mechanism to evade the most stringent experimental constraints on the strength of the field's coupling to matter. For 
example: in the solar system, the chameleon can be very light thus mediate a long-range force. The limits on such 
forces are very tight, [IEUt]]- However, since the chameleon only couples to a small fraction of the matter in large 
bodies i.e. that fraction in the thin-shell, the chameleon force between the Sun and the planets is very weak. As a 
result the chameleon has no great effect on planetary orbits, and the otherwise tight limits on such a long-range force 
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Sketch of chameleon mechanism: Low Density Background 



Sketch of chameleon mechanism: High Density Background 
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FIG. 1: Sketch of the chameleon mechanism for a runaway potential: V ~ </> -4 - The sketch on the left is for a low density 
background, whereas the drawing of the right shows what occurs when there is a high density of matter in the surroundings. We 
can clearly see that the position of the effective minimum, 4> c , and the steepness of the effective potential near that minimum, 
depends on the density. A shallow minimum corresponds to a low chameleon mass. The mass of the chameleon can be clearly 
seen to grow with the background density of matter. 



Sketch of chameleon mechanism: Low Density Background 



Sketch of chameleon mechanism: High Density Background 





FIG. 2: Sketch of the chameleon mechanism for a potential with a minimum at (j> — 0: V ~ <f> . The sketch on the left is for 
a low density background, whereas the drawing of the right shows what occurs when there is a high density of matter in the 
surroundings. We can clearly see that the position of the effective minimum, <j> c , and the steepness of the effective potential 
near that minimum, depends on the density. A shallow minimum corresponds to a low chameleon mass. The mass of the 
chameleon near 4> c can be clearly seen to grow with the background density of matter. 

are evaded, [2||. In section UTT1 we will show that the presence of a thin-shell effect is intimately linked to non-linear 
nature of chameleon field theories. 



B. Chameleon to Matter Coupling 



When a scalar field, </>, couples to a species of matter, the effect of that coupling is to make the mass, m, of that 
species of particles ^-dependent. This can happen either at the classical level (i.e. in the Lagrangian) or a result of 
quantum corrections. We parameterise the dependence of m on by 



m(4>) = toqC 



M p i 



(1) 
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where M p i is the Planck mass and mo is just some constant with units of mass whose definition will depend on one's 

choice of the function C (^Jj ■ P defines the strength of the coupling of 4> to matter. We shall say more about the 

definition of (3 below. A (^-dependent mass will cause the rest-mass density of this particle species to be (f> dependent, 
specifically 

ptf)=p c(£P\. (2) 



The coupling of 4> to the local energy density of this particle species is given by: dp{<j))/d(j) which is: 

dptf) = pi ( ^ 
d</> a \M pl J M pl ' [6> 

where B(x) = \nC(x) and B'(x) — dB{x) / dx. Throughout this work we will, for simplicity, assume that our 
chameleon field, (j>, couples to all species of matter in the same way, however we will keep in mind the fact that, 
generically, different species of matter will interact with the chameleon in different ways. We shall see in section [Villi 
that, if C, and hence B, are at least approximately the same for all particle species, then cosmological bounds on 
chameleon theories will require that \/3B'(/3(f>/Mpi)(f>/M p i\ < 0.1 everywhere since the epoch of nucleosynthesis. We 
preempt this requirement and use it to justify the linearisation of B{j3<f>/M p i): 

For this to be a valid truncation we require (B"(0)/B'(0))f34>/M pl < 1. So long as \B"(0)\ < 10|S'(0)|, the cosmolog- 
ical bounds on (j) will then ensure that the above truncation of the expansion of B is a valid one. The only forms of B 
that are excluded from this analysis are the ones where |B"(0)| > 10|B'(0)|; we generally expect B"(0) ~ O(B'(0)). 

Provided B'(0) ^ 0, we can use the freedom in the definition of (3 to set B'(0) = 1. When this is done, (3 quantifies 
the strength of the chameleon-to-matter coupling. 

For example: a particular choice for C that has had some favour in the literature, [13, El, is C = e^^ 1 ? 1 for some 
k. It follows that B = k(f>/M p i, and so we choose (3 = k which ensures -B'(O) = 1, B"(0) = 0. 

We wish to construct our chameleon theories to be compatible with Einstein's conception of gravity. By this we 
mean that we wish them to display diffeomorphism and Poincare invariance at the level of the Lagrangian. A natural 
consequence of Poincare invariance is that the chameleon couples to matter in a Lorentz invariant fashion. For a 
perfect fluid this implies that cf> will generally couple to some linear combination of p and the fluid pressure (P) i.e 
P + luP. The simplest way for the chameleon to couple to matter in a relativistically invariant fashion is for it to couple 
to the trace of the energy momentum tensor; in this case u> = —3. This said, apart from in the early universe and in 
very high density objects such as neutron stars, P/p <C 1, and so the precise value of w is not of great importance. 
Apart from where such an assumption would be invalid, we will take P/p <C 1 and set P = 0. 



C. A Lagrangian for Chameleon Theories 

It is possible to couple the chameleon to matter in a number of different ways, and as such it is possible to construct 
many different actions for chameleon theories. A reasonably general example of how the chameleon can couple to 
trace of the energy momentum tensor is given by following Lagrangian density 

+ £ m (^\ 5 «), (5) 

where C m is the Lagrangian density for normal matter. This Lagrangian was first proposed in ref. [25j ] . The index i 
labels the different matter fields, ip^\ and their chameleon coupling. The metrics g$, are conformally related to the 
Einstein frame metric g^ u by g[ll = Q^g^ where = C \i)((3^ '4 1 1 M p i) . The C^(-) are model dependent functions 
of f3^<t)/M p i. The /3« are chosen so that B'^(0) := (lnC(j))'(0) = 1. R(g) is the Ricci-scalar associated with the 

Einstein frame metric. For simplicity we will restrict ourselves to a universal matter coupling i.e g^l — g tlVl Gu\ = C 
and /3® = (3. We define T^ v = (2/ y /^g)5£ m /5g lll/ . It follows that T = T^g^ = p - 3P, where p is the physical 
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energy density and P is the sum of the principal pressures. In general p and P are ^-dependent. With respect to this 
action, the chameleon field equation is 

-U<t> = V <t ,{4>) + -f . (6) 

As mentioned above, \(3(f>/M p i\ < 0.1 is required for the theory to be viable and so it is acceptable to approximate 
B'(j3(j)/Mpi) by B'(0). We then scale (3 so that B'(0) = 1. The requirement, \j3(j>/M p i\ < 0.1, also ensures that p and 
P are independent of <f> to leading-order. The field equation for <f> is therefore 

-«. W1 + ffl»^. (7) 

Mpi 

The above Lagrangian should not be viewed as specifying the only way in which <f> can couple to matter. When 
one considers varying constant theories, the matter coupling often results from quantum loop effects [95|. However, 
despite the fact that many different Lagrangians are possible, it is almost always the case that the field equation for 
<f> takes a form very similar to the one given above. 



D. Intrinsic Chameleon Mass Scale 



(3 quantifies the strength of the chameleon coupling. := M p i/ (3 can then be viewed as the intrinsic mass scale of 
the chameleon. Although precise calculations of scattering amplitudes fall outside the scope of this work, we expect 
that chameleon particles would be produced in large numbers in particle colliders that operate at energies of the order 
of or greater. 

It is generally seen as 'natural', from the point of view of string theory, to have M p i. When this happens 

the chameleon has the same energy-scale as gravity. It has also been suggested that the chameleon field arises from 
the compactification of extra dimensions, [1], if this is the case then there is no particular reason why the true Planck 
scale (i.e. of the whole of space time including the extra-dimensions) should be the same as the effective 4-dimensional 
Planck scale defined by M p i. Indeed having the true Planck scale being much lower than M p i has been put forward as 
a means by which to solve the Hierarchy problem (e.g. the ADD scenario [MEM Hoi). In string-theory too, there 
is no particularly reason why the string-scale should be the same as the effective four-dimensional Planck scale. It is 
also possible that the chameleon might arise as a result of new physics with an associated energy scale greater than the 
electroweak scale but much less than M p i. In light of these considerations it would be pleasant if = M p i/j3 <C M p i, 
say of the GUT scale, or, if we hoped to find traces of it at the LHC, maybe even the TeV-scale. 

A positive detection of a chameleon field with such a sub-Planckian energy scale could provide us with the first 
evidence for new physics beyond the standard model, but below the Planck scale. 

As pleasant as it might be to have M$ <C M p i, it is generally agreed that the current experimental bounds on 
the existence of light scalar fields rule out this possibility fl3l . [lil Il5l Indeed, in the absence of a chameleon 



mechanism similar to that proposed in [18 
coming from Lunar Laser Ranging (LLR), 



, 25], bounds on the violation of the weak equivalence principle (WEP) 
Mill, limit \/3\ < 10~ 5 for a light scalar field. This implies Mj, > M pl . 



If the Planck scale is supposed to be associated with some fundamental maximum energy, such a large value of 
seems highly unlikely. Even if a (non-chameleon) scalar field has a mass of the order of 1mm -1 , then we must have 
(3 < 10-\ [13. 

One of the major successes of the proposal of chameleon field by Khoury and Weltman, [l?], Hf|, was that chameleon 
fields can, by attaining a large mass in high density environments such as the Earth, Sun and Moon, evade the 
experimental limits coming from LLR and other laboratory tests of gravity. In this way, it has been shown the scalar 
fields in theories that possess a chameleon mechanism can couple to matter with the strength of gravity, (3 ~ O(l) 
and still coexist with the best experimental data currently available. Even though (3 ~ 0(1) has been shown to be 
possible, (3 3> 1 is still generally assumed to be ruled out. 

In this paper, however, we challenge this assumption and show that it is indeed feasible for (3 to be very large. 
Moreover » Mqut ~ 10 15 GeV or ~ 1 TeV are allowed. Tantalisingly the experimental precision required to 
detect such a sub-Planckian chameleon theory is already within reach. Large matter couplings are allowed in normal 
scalar field theories but only if the scalar field has a mass greater than (0.1 mm) -1 . This is not the case for chameleon 
theories. We shall show that the mass of the chameleon in the cosmos, or the solar system, can be, and generally is, 
much less than 1 m — . 



8 



E. Initial Conditions 



Even though the term chameleon field sounds rather exotic, in a general scalar field theory with a matter coupling 
and arbitrary self-interaction potential, there will generically be some values of <f> about which the field theory exhibits 
a chameleon mechanism. Whether or not <fr ends up in such a region will depend on its cosmological evolution and 
one's choice of initial conditions. The importance of initial conditions was discussed in [l8[. In that paper the potential 
was chosen to be of runaway form V oc 4>~ n , n > 0. We will review what is required of the initial conditions for such 
potentials in section IVIIII below. We shall see that the larger f3 is, the less important the initial conditions become. 
We will also see that the stronger the coupling, the stronger the chameleon mechanism and so the more likely it 
becomes that a given scalar field theory will display chameleon- like behaviour. This is one of the reasons for wanting 
to have a large value of (3. 



F. The Importance of Non-Linearities 



Chameleon field theories necessarily involve highly non-linear self-interaction potentials for the chameleon. These 
non-linearities make analytical solution of the field equations much more difficult, particularly when the background 
matter density is highly inhomogeneous. Most commentators therefore linearise the equations of chameleon theories 
when studying their behaviour in inhomogeneous backgrounds 0, HO, HH, H2] • Such an approximation may mislead 
theoretical investigations and result in erroneous conclusions about experiments which probe fifth force effects. In this 
paper we shall show, in detail, that this linearisation procedure is indeed very often invalid. When the non-linearities 
are properly accounted for, we will see that the chameleon mechanism becomes much stronger. It is this strengthening 
of the chameleon mechanism that opens up the possibility of the existence of light cosmological scalars that couple 
to matter much more strongly than gravity (/3 3> 1). 



G. The Chameleon Potential 



The key ingredient of a chameleon field theory, in addition to the chameleon-to-matter coupling, is a non-linear 
and non-quadratic self-interaction potential V{4>). It has been noted previously that V{<j>) could play the role of an 
effective cosmological constant [18j |. There are obviously many choices one could make for V(<p), and whilst we wish 
to remain suitably general in our study, we must go some way to specifying V((f>) if we are to make progress. On e 
quite general form that has been widely used in the literature is the Ratra-Peebles potential, V{4>) = M (M/<f) n 
where M is some mass scale and n > 0; chameleon fields have also been studied in the context of V(<j)) = kef) 4 /4! 
In this paper we will consider both of these types and generalise a little further. We take: 

V{cf>) = AAf 4 (M/0)", (8) 

where n can be positive or negative and A > 0. If n ^ —4 then we can scale M so that without loss of generality 
A = 1. When n — —4, M drops out and we have a </> 4 theory. When n > this is just the Ratra-Pebbles potential. 



34]. 
35|. 



H. Chameleon Field Equation 



With these assumptions and requirements, the chameleon field, </>, obeys the following conservation equation: 

For this to be a chameleon field we need the potential gradient term, = —n\M 3 (M / <fi) n+1 , and the matter coupling 
term, /3(p + u>P)/M p i, to be of opposite signs. It is usually the case that j3 > and P/p <C 1. If n > we must 
therefore have (f> > 0. In theories with n < we must have cj> < and n — —2p where p is a positive integer. We must 
also require that the effective mass-squared of the chameleon field, = V^, be positive, non-zero and depend on 
4>. These conditions mean that we must exclude the region —2 > n > 0. If n = —2, n = —I or n = then the field 
equations for <fi would be linear. 
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I. Natural Values of M and A 

When n ^ —4, one might imagine that our choice of potential has arisen out of an expansion, for small (M /4>) n , of 
another potential W(<fi) — M 4 f ((M / 4>) n ) where / is some function. We could then write: 

W ~ Af 4 /(0) + M\f'(0) (^j , (10) 

where M is some mass-scale. We define M so that the second term on the right hand side of the above expression 
reads M 4 (M/<fi) n . The first term on the right hand then plays the role of a cosmological constant M 4 /(0) « p\- 
Assuming that both /(0) and f'(0) are 0(1), we would then have M « M w (oa) 1//4 « (0.1mm) -1 . It is for this 
reason that one will often find (0.1 mm) -1 referred to as a 'natural' value for M, [3 EE]- When n = —4 we naturally 
expect A w 1/4!, [H]. 



III. ONE BODY PROBLEM 



In this section, we consider the perturbation to the chameleon field generated by a single body embedded in 
background of uniform density pb- For simplicity we shall model the body to be both spherical and of uniform density 
p c . This analysis will prove vital when we come to calculate the force between two bodies that is induced by the 
chameleon field. We assume space-time to be Minkowski (at least to leading order) and we also assume that everything 
is static. Under these assumptions □ — > —V 2 where V 2 </> = r~ 2 (r 2 </>')'; <f>' = d(j>/dr. Whilst this problem has been 
considered elsewhere in the literature [13, EH HH , most commentators have chosen to linearise the chameleon field 
equation, eq. ©, before solving it. This linearisation is, however, often invalid. In this section, we begin by briefly 
reviewing what occurs when it is appropriate to linearise eq. ([9]), and, in doing so, note where the linear approximation 
breaks down. In some cases, even though it is not possible to construct a linearised theory that is valid everywhere, 
we shall demonstrate, using the method of matched asymptotic expansions, how to construct multiple linearisations of 
the field equation, each valid in a different region, and then match them together to find an asymptotic approximation 
to <j) that is valid everywhere. When this is possible, the chameleon field, cf>, will behave as if these were the solution 
to a consistent, everywhere valid, linearisation of the field equations; for this reason we deem this method of finding 
solutions to be the pseudo-linear approximation. If a body is large enough, however, both the linear and pseudo-linear 
approximations will fail. We shall see that, when this happens, <p behaves in a truly non-linear fashion near the surface 
of the body. The onset of non-linear behaviour is related to the emergence of a thin-shell in the body. The linear 
approximation is discussed in section fill Al whilst the pseudo-linear approximation is considered in section fill Bl Wc 
discuss the non-linear regime in section fill CI 

We take the body that we are considering to be spherical with radius R and uniform density p c . Assuming spherical 
symmetry, inside the body (r < R), (j) obeys: 

2d0 ^n f3 f M Y +1 , PPc 



and outside the body, (r > R), we have: 



d2 ^ 2 - d 4 = -nXM*( A 4Y +1 + ^. (12) 



dr 2 r dr \ <fi J ' M p i 

The right hand side of eq. (Ilip vanishes when <f> = <f> c where 

PP. 



= M 



iXMpiM 3 



This value of cj> corresponds to the minimum, of the effective potential of the chameleon field, inside the body. Similarly, 
the right hand side of eq. (|12p vanishes when <fi = <fi b where 



M 



nXMpiM 3 ] 
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This value of <fi corresponds to the minimum, of the effective potential of the chameleon field, outside the body. For 
large r we must have tfi w (fib- Associated with every value of <fi is an effective chameleon mass, ((f)), which is the 
mass of small perturbations about that value of 4>. This effective mass is given by: 

/ /I f \ "+ 2 

™W = VfflM = »(" + 1)AM 2 . (13) 

We define m c = m^(cfi c ) and mj, = m^, ((/>&). We shall see below that the larger the quantity m c R, the more likely it is 
that a body will have a thin-shell. In this section we shall see both why this is so, and precisely how large m c R has 
to be for a thin-shell to appear. Throughout this section we will require, as boundary conditions, that 

= 0. 



dr 



= and — — 

dr 



A. Linear Regime 



We assume that it is a valid approximation to linearise the equations of motion for <fi about the value of <f> in the far 
background, (fib- For this to be possible we must require that certain conditions, which we state below, hold. Writing 
<fi = 4>b + <Pi i t ne linearised field equations are: 

-l 



d 2 (t>i 
dr 2 



2d(j) 1 
r dr 



-nM i 



M 



P(Pc - Pb) 



M, 



pi 



H(R 



(14) 



where H(R — r) is the Heaviside function: H(x) = 1, x > 0, and H(x) — 0, x < 0. 
potential to be valid we need: 



(3p b 
M pl ' 

For this linearisation of the 



V, 



< 1. 



This translates to \(fix/(fib\ < \n 
implies that: 



. Also, for this linearisation to remain valid as r 



oo, we need <fii — > 0, which 



n+l 



PPb 



Defining Ap c — p c — p^, and solving the field equations, we find that outside the body (r > R) we have: 

PAp c e mb{R - r) /tanh(m 6 i?) -m b R" 



Mpim 



nibr 



\ 1 + tanh(m6-R) 



Inside the body (r < R), (f>i is given by 



(3Ap c f3Ap c (1 + m b R) e- m » H sinh(m 6 r) 



m%M p i M v im% 



m,br 



The largest value of \$\l<$>b\ occurs at r = and so, for this linear approximation to be valid, we need: \4>i(r 
0)/(j>b\ < \n+ 1|~ . This requirement is equivalent to the statement that 



1(1 + rribR)e 



i| ^fi ^ - m 2 R 2 Ht£i = ^fi ( Eh 



Pb 



Pb 



^ (m c R) 2 



< 1. 



where '~' means "asymptotically in the limit rribR — * 0". It is often the case that the background is much less dense 
than the body i.e. pb <C p c . If this is the case then it is clear, from the above expression, that there will be a distinct 
difference between theories with n > and those with n < —4. In theories with n > 0, the lower the density of the 
background, the better the linear theory approximation will hold, whereas when n < —4 the opposite is true. This 
can be understood by considering the relation (fib oc p b 1 '^ n+1 ', If n > 0, the smaller pb becomes, the larger (fib will 
be. It is therefore possible for larger perturbations in <fi to be treated consistently in terms of the linearised theory. 
If, however, we have that n < — 4 then (fib — * as pb — > and the opposite is true. 

We can, however, use the method of matched asymptotic expansions to show that the region where behaviour, 
similar to that which would be predicted by linearised theory, occur is significantly larger than one would have 
guessed simply by requiring that the linear approximation hold. 

The results of this section, as well as those of sections IIII Bl and IIII C[ are summarised in section IIII Dl below. 
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B. Pseudo-Linear Regime 



The defining approximation of the pseudo-linear regime (for both positive and negative n) is that inside the body: 

n+l 

< 1 



(r) 

This is equivalent to f3p c /M p i 3> nM 3 (M / fo[r)) n+1 . When this holds we find 4> ~ 4>( r ) inside the body, where this 
defines for) and: 

1 d 2 (rfo = fa 
r dr 2 M p i ' 

It follows that 

Pp c r 2 



! >o 



6M p i 



In this case l ~' means "asymptotically as (fo/ for)) n+1 — > 0". Outside the body, we can find a similar asymptotic 
approximation: 



" 2M p i 3M pl r 

For this to be valid we must ensure that the neglected terms, in the above approximation to fo are small compared 
to the included ones; this requires that: 



3 Jo Jo \<Kr"). 
For large r we expect, as we did in the previous section, that <f> —> fo, and so: 

fo~ fo* =fo , 

r 

which will remain valid whenever Ae~ mbT 'for <C 1 1 / (n, + 1). We shall refer to <f> as the inner approximation to <f>. 
Similarly, fo is the outer approximation. So far both A, and the value of fo, remain unknown constants of integration. 
In general, when <j) ~ fo" we will not also have <f> ~ <j) (and vice versa). If, however, there is some intermediate region 
where both the inner and outer approximations are simultaneously valid, then we can match both expressions in that 
intermediate region and determine both fo and A, [43l |44| . A detailed explanation of the use of matched asymptotic 
expansions with respect to cosmological scalar fields is given in [45j . 

For the moment we shall assume that such an intermediate region does exist. We check what is required for this 
assumption to hold in appendix [X] and present the results of that analysis below. Given an intermediate region, we 
find: 

A = P pcB? 
3M pl ' 

, PPcR 2 Pp c m b R 3 
2M p i ZM p i 

The external field produced by a single body in the pseudo-linear approximation is: 

(3p c R 3 e' mbr fo {m c R) 2 e- mbr 

~ 0b 3M pl r * fo~ fo 3(n+l)(r/Ry 

and the field inside the body is given by: 

(3p c R 2 (3p c m b R 3 (3p c r 2 

2M„; 3M„i 6M„; 
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In appendix ^ we find that for the pseudo- linear approximation to hold, we must require 
m c R « rmn ((18)V6(^)^ ) V2|^TlT|(^)" 

m c R « mm(V3(^)^,^^|fe)^ 
m c i? < 1, n = -4. 

When n — —4, we actually find a slightly different asymptotic behaviour of <j) outside the body, precisely: 



1 



, n < -4, 
n > 0, 



1 



2A(y + ln(min(r/i?, l/m&))) 2r ' 



(17a) 

(17b) 
(17c) 



(18) 



where yo is an integration constant and: 



(m c R) 3 



2yo 



3/2 



The conditions given by eqs. (|17b-c) ensure that the pseudo-linear approximation is everywhere valid. When these 
conditions fail, non-linear effects begin to become important near the centre of the body. As m c R is increased further 
the region where non-linear effects play a role moves out from the centre of the body. Eventually, for large enough 
m c R, the non- linear nature of chameleon potential, V(4>) 7 is only important in a thin region near the surface of the 
body; this is the thin-shell. Since the emergence of such a thin-shell is linked to non-linear effects becoming important 
near the surface of the body, it must be the case that the assumption that (f> is given by equation (fl~5f (or by eq. (fl8|) 
in n = —4 theories) breaks down for some r > R. By this logic, we find, in appendix \X[ that a thin-shell occurs when: 



m c R > min^lS) 1 / 6 ^)™,^ 

^3 3< " +2) , VW^+M 

m r R > 4, n = -4. 



n < 



n > 0, 



(19a) 

(19b) 
(19c) 



In both eqs. (fTTb-b) and (fT9"k -b) the second term in the min( • , • ) is almost always smaller than the first when 
Pb/Pc < 1 ^ mb/m c <C 1. The behaviour of <p both near to, and far away from, a body with thin-shell is discussed 
below in section UlI CI The results of this section are summarised in section Ull Dl Note that the thin-shell conditions 
, eqs. (fT9"k-c), necessarily imply that m c R ^> 1. 



C. Non-linear Regime 



We have just seen, in eqs. (flUb-c). that for non-linear effects to be important, and the pseudo- linear approximation 
to fail, we must have m c R ^> 1. In this regime the body is, necessarily, very large compared to the length scale l/m c . 
We expect that all perturbations in (f> will die off exponentially quickly over a distance of about 1 / m c and, as such, 
(j) ps (j) c will be almost constant inside the body. Any variation in the chameleon field, that does take place, will occur 
in a 'thin-shcir of thickness AR w l/m c near the surface of the body. It is clear that m c R ^> 1 implies AR/R <C 1. 
In this section we will consider both the behaviour of the field close to the surface of the body, and far from the body. 



1. Close to the body 

We noted above that m c R ^> 1 implies AR/R C 1, we shall demonstrate this is a rigorous fashion below. Given 
AR/ R <C 1 , when we consider the evolution of </> in the thin-shell region, we can ignore the curvature of the surface 
of the body, to a good approximation. 

We therefore treat the surface of the body as being flat, with outward normal in the direction of the positive ir-axis. 
The surface of the body defined to be at x = (i.e. x = r — R). Since the shell is thin compared to the scale of the 
body, we are interested in physics that occurs over length-scales that are very small compared to the size of the body. 
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We therefore make the approximation that the body extends to infinity along the y and z axes and also along the 
negative x axis. Given these assumptions, we have that <f> evolves according to 



= -n\M 3 (M/cf>) n+1 +/3 Pc /M pl 

As a boundary conditions (BCs) we have <j> — ► <p c and d<fi/dx — > as x — ► • 
of the above equation is: 



-co. With these BCs, the first integral 



1 fd(j> 

2 V dx 



Outside of the body, we assume that 
that: 



AM 4 [(M/0)" ~ (M/4, c ) n ] + f^(0 ~ 0c). (20) 

M p i 

-> <j)b as x — > oo, and that the background has density pb p c - Assuming 



dx 2 



» 



2dcj) 
r dx 



then we can ignore the curvature of the surface of the body and, in x > 0, we have: 

Uf) 2 = \M\M/W ~ AM 4 (M/0 b )» + %?L" 



2 V dx 



pi 



(21) 



where 0{, is as we have defined above. Our assumption that \d 2 (f>/dx 2 \ 3> (2/r)d0/dx then requires that: 



2V2 



m^(0)r V n 



n + 1 



l-(n+l)(£ 



<Pi, 



n+1 



1/2 



< 1. 



Provided the pseudo-linear approximation breaks down, and that the body has a thin-shell, we expect that, near the 
surface of the body, ~ 0(0 C ). It follows that, whenever p c ^> pb, (0/0b)™ -C 1 and (0/0b)™ +1 <C 1. The above 



condition will therefore be satisfied provided that m c R y/8(n + l)/n; this is generally a weaker condition than the 
requirement that the body satisfy the thin-shell conditions, eqs (|T5k -c). On the surface at x = 0, both and d0/dx 
must be continuous. By comparing the expressions for d0/dx inside and outside the body we have: 



0(0) 



1 

n 



We now check that we do indeed have a thin shell i.e. AR <§; R. We expect that, near the surface of the body, 
almost all variation in will concentrated into a shell of thickness AR. We define m sur f by: 

d4> 



^ X ~dx = ^ X = °) ~ ^ c = ^ c / n = msurf ~dx~( X = 
m mrf 1S then, approximately, the length scale over which any variation in dies off. As it happens, m sur f is also 



«i o r„f- For this shell to be thin, and for us to be 

sal j 



the mass of the chameleon field at x = 0. It follows that AR ~- ' 

justified in ignoring the curvature of the surface of the body, we need AR/R <C 1 or equivalently m sur fR 3> 1. We 
find (assuming pb <C p c ) that: 



TUsur f R ' 



n + 1 



n/2+1 



m c R~ 0(m c R), 



and so m sur fR 3> 1 follows from m c R 3> 1, and AR ~ ©(rn^ 1 ). m sur fR 3> 1 will be automatically satisfied whenever 
the thin-shell conditions eqs. (fT9a -c) hold. 

Whenever p c ^$> pb, eq. [2TI will . near r = R, be well-approximated by: 



1 /d0 

2 \ dx 



AM 4 (M/0)". 
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Solving this under the boundary conditions <f>(x = 0) = (1 + \/n)<j) c and 4>/4>c — * 4>b/4>c ~ as x — > oo we find 

1 „ |n + 2|(r-.R) | /n + l\ " /2+1 J_ 
m^,(</>) v /2n(n + 1) V " / TO c 

This approximation will therefore break-down when m^fyr ~ 0(1), which occurs when r — i? ~ O(R). We can see 
that, if r - i? > ^/2n(n + l)/(|n + 2|m s ur f) then m^, and hence also <p, will be independent of m c and hence also of 
<p c and /3 at leading order. Since m sur fR 3> 1, there will be some region where eq. (|22l) is both valid and, to leading 
order, independent of (3. 

Although, in this approximation, we cannot talk about what occurs for (r — R) > R, it seems likely, in light of 
the behaviour seen when (r — R) <C R, that, whenever r 3> l/ m «w/ m l/m c , the perturbation in 0, induced by an 
isolated body with thin-shell, will be independent of the matter coupling (3. We confirm this expectation in section 
IIIIC 21 below. 



2. Far field of body with thin-shell 



We found above that the emergence of a thin-shell was related to non-linear effects being non-negligible near the 
surface of the body. We noted that a thin-shell will exist whenever conditions (|19a -c) hold. However, even when 
these conditions hold, we do not expect non-linear effects to be important far from the surface of the body. Indeed, 
for large r we should expect that <f> takes a functional form similar to that found in the pseudo-linear approximation 
i.e. as given by eq. (JTSJ) (or eq. {T5J) if n = —4). Although the functional form should be similar, in order to find 
the correct behaviour, one must replace (m c R) in equations (fT5|) and (fTgj) by some other quantity C, say. We show 
that, to leading order as r — > oo, C is independent of the matter coupling (3 and the density of the central body p. 
This confirms the expectation of section IIII C II above. The analysis for n = —4 is slightly more involved than it is 
for other values of n. We therefore consider the n = —4 case separately below and in appendix [S] The analysis 
for theories with runaway potentials that become singular at <p = (i.e n > theories) is much simpler than it is 
for theories where the potential has a minimum at <j> = and which are non-singular for all finite <f> i.e. (n < —4 
theories): we therefore consider the n < and n > cases separately. 



Runaway Potentials (n > 0) 

Away from the surface of the body we expect that non-linear effects will be negligible and as r — > oo we will have: 

£) e -m b r 



for some D where <j>b and mb are the values of the chameleon and its mass in background. It is clear from the field 
equations however that V 2 < V 2 0^°- ) and so, given the boundary conditions at infinity, (f> < </>(°) outside the body. 
In n > theories there is a singularity of the potential, and hence also of the field equations, at <j> = 0. It is clear that 
this singularity cannot be reached in any physically acceptable evolution and so we must always have > 0, which in 



turn implies (f)^ ' > outside the body. The minimum value of 
have: 



outside the body occurs at r = R and so we must 



D < <f> b e mbR R. 

In most cases of interest iribR <C 1 and so we have: 

D < (f> b R. 

This upper bound on D defines a critical form for the field outside the body: 

e m b (R-r) ft> 



1 - 



No matter what occurs inside the body (r < R) we must have <f> > 4> C rit outside the body as r — > oo. This implies 
that: 



dr 



(1 + m b r)(t)be mb{ - R - r)R 4> b R 
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as r — > oo. Ignoring non-linear effects, 4> > 4> cr it is satisfied by all bodies that satisfy the conditions for the pseudo- 
linear approximation (eqs. I|17b)) but would be violated, in the absence of non-linear effects, by those that satisfy the 
thin-shell conditions (eqs. (fT9|) ). We must therefore conclude that non- linear effects near the surface of body with 
thin-shells ensure <p > <j) cr it is always satisfied as r — > oo. Furthermore, if </> 3> 4> C rit then 

d4> 
dr 

and it follows from section fill B I that the pseudo-linear approximation is valid for all r, which further implies that the 
body cannot have a thin-shell. Thus thin-shelled bodies must actually have (f> being only greater less than (frcrit as 
r — > oo. We are therefore justified in using 4> cr it to approximate the far field of a body with a thin-shell. In summary: 
in n > theories, the far field of a body with a thin-shell has the following form: 

<\> ~ fa - fae m ^ R ~ r)R r. 

We note that this form, and the arguments with which we have derived it, do not depend, in any way, on the physics 
inside r < R. The critical form of <p is determined entirely by the form of the potential and the background value of 

<$>■ 



kR 



Potentials with minimum (n < 0) 

For n > theories the singularity of the potential at <fi = allowed us to determine asymptotic form of tf> outside 
a body with a thin-shell. In n < theories, however, the potential is well-defined for all finite <fi and so we cannot 
play the same trick as we did above. The n = —4 case is special and treated in great detail in appendix [5] When 
n = —4, we find that the far field of a body with thin-shell is given by: 



ry/2\(l + 41n(min(r/i?, l/m b R))) 2ry/2X ln(min(r/E, l/m b R)) 

A thin-shell is certainly present whenever m c R > 4. For other negative values of n we use a semi-analytical method. 
We saw when deriving the thin-shell conditions for n < theories that the background value of <j> plays only a 
negligible role since 4>/fa ^ 1 near the body and, in most cases, tyi^R <ti 1. Assuming m b R <C 1, we simplify our 
analysis by setting </>;, = 0. Far from the body non- linear effects are sub- leading order and we expect: 

= + o(l/r). 

r 



We now define a new coordinate s = \J\n\A~t n+2 } Mr and u = —<f>/AM for some constant A. With these definitions 
the full field equation for 4> outside the body (with fa = 0) becomes: 

1 d ( s2 du\ =u _ n _ 



s 2 ds \ ds J 
and as s — > oo: 

DA-^ 

U —— + o(l/8). 

We set A n+i = D 2 /\n\ so that u~l/s and define t = 1/s so that the field equations become: 

d 2 u _ u-"- 1 

~d¥ ~ [26) 

The asymptotic form of u as r — » oo, t — > 0, requires that u(t = 0) = and du/dt(t = 0) = 1 exactly. With these 
boundary conditions we numerically evolve eqn. ()23|1 towards larger t (smaller r). As one might expect from such 
an elliptic equation, with these boundary conditions, a singularity occurs at some finite t which we label t ma x- We 
use our numerical evolutions to determine t max for each n. For the evolution of <p to remain non-singular up to the 
surface of the body, we need r = R to correspond to a value of t < t max . The limiting case is given by t(R) — t max . 
This limiting case determines a critical form for the (f> field which occurs when A = A min where: 



nt 2 max \ — (MR)- 
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TABLE I: Values of y(n) = t 



n 


7(n) 


-12 


14.687 


-10 


10.726 


-8 


6.803 


-6 


3.000 




This corresponds to the following critical asymptotic form for <p: 

„ +4 e ~ m br 

(MR)*+s , 

r 

where we have reinserted the (almost always negligible) 4>b and m& dependence. We use numerical integration 
to calculate the value of 7(n) := tmax for different values of n. Our results are displayed in table Q] Physically 
acceptable non-singular evolution implies that asymptotically 4>/4> C rit < 1- If ffr/ficrit <C 1 then the conditions 
of the pseudo-linear approximation are satisfied and so the body cannot have a thin-shell. Thin-shelled bodies 
must therefore almost saturate this bound on </> and so <fi cr it provides a good approximation to the asymptotic 
behaviour cf> outside thin-shelled bodies. We note that, as in the n > case, the existence of a critical form for </> de- 
pends in no way on the what occurs inside the body and, as such, is independent of both /3 and the density of the body. 

Critical Behaviour 

The existent of a critical form for <p when r » i! implies that, no matter how massive our central body, and no 
matter how strongly it couples to the chameleon, the perturbation it produces in tf> for r>]J takes a universal value 
whenever the thin-shell conditions, eqns (fT§k-ch hold. 

When n ^ —4, the critical form of the far field, depends only on M, n, R and on the chameleon mass in the 
background, mf,. When n = —4 the critical form for the far field depends only on A, R and m^. For all n, the far 
field is, crucially, found to be independent of the coupling, /3, of the chameleon to the isolated body. This is one 
of the main reasons why /3 3> 1 is not ruled out by current experiments. The larger (3 becomes, the stronger the 
chameleon mechanism and so the easier it is for a given body to have a thin-shell. However, the far field of a body 
with a thin-shell is independent of f3, and so, in stark contrast to what occurs for linear theories, larger values of f3 
do not result in larger forces between distant bodies. Defining the mass of our central body to be M. = 47r / o c i? 3 /3 we 
can express this critical behaviour of the far field in terms of an effective coupling, /3 e //, defined by: 

, , f3 eff Me- m > r 

4> ~ 4>b - 



AttM p1 7 

when r R. Assuming Pb/p c 1 we find that: 



M \ \n 

Peff = ~M MR { ml J ' n> ° (24b) 
/3 eff (r) = ^^(2Xln(mm(r/R,l/m b R))y 1/2 , n = -4. (24c) 

The j3 independence of (3 e ff was first noted, in the context of <p A theory, in [37| . However, the authors were mostly 
concerned with region of parameter space f3 < 1, A <^ 1; in our analysis we go further: considering a wider range 
of theories and also the possibility that (3 3> 1. /3-independence was also present in the original work of Khoury 
and Weltman [13, for n > theories. However, in those works, the (3 independence together with its important 
implications for experiments, was not commented on. Especially those that search for WEP violations were not 
considered. As we shall see in section IVII below, this (3 independence means that if one uses test-bodies with the 
same mass and outer dimensions then in chameleon theories, no matter how much the weak equivalence principle is 
violated at a particle level, there will be no violations of WEP far from the body. Simply because the far field is 
totally independence of both the body's chameleon coupling and its density. 
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In this work, we have also shown that this (3 independence is a generic feature of all V oc (p~ n chameleon theories 
and it is not simply as artifact of the runaway (n > 0) potentials considered in [T3, H^|. Indeed there are good 
reasons to believe that similar behaviour will be seen in chameleon theories with other potentials. As we mentioned 
in the introduction, the field equations for chameleon theories are necessarily non-linear. It is well-known that, that 
in non-linear theories, the evolution of arbitrary initial conditions will gencrically be singular. If one wishes to avoid 
singularities then tight constraints on the initial conditions must be satisfied. When considering the field outside an 
isolated body, these conditions will, generally, require that |d</>/dr| is smaller than some critical, r-dependent, value. 
As a result, there will a critical, or maximal, form that the field produced by a body can take. This precisely what we 
have found for <\T n theories. The form of this critical far field will depend on the nature of the non-linear potential, 
and possibly the coupling of </> to any background matter, but, since we are outside the body, it cannot depend on the 
coupling of the chameleon to the body itself. Again, this is precisely what we have seen for </>~™ chameleon theories. 

We can understand the /3-independence, in a slightly different way, as follows: just outside a thin-shelled body, the 
potential term in eq. ([9]) is large and negative (~ 0(—(3p/M p i)), and it causes <p to decay very quickly. At some point, 
<j) will reach a critical value, (f> cr it, that is small enough so that non-linearities are no longer important. Since this all 
occurs outside the body, <j) cr it can only depend on the size of the body, the choice of potential (M, A, n) and the mass 
of 4> m the background, mi,. This is precisely what we have found above. 

We have seen above that the far field of a body with thin-shell is independent of the microscopic chameleon-to-matter 
coupling, (3. This is one of the vital features that allows theories with (3 ^> 1 to coexist with the current experimental 
bounds. It is also of great importance when testing for WEP violations, since any microscopic composition dependence 
in (3 will be invisible in the far field of the body. We discuss these issues further in section fVTl where we consider the 
experimental constraints on f3, M and A in more detail. 

The results of this section are summarised in section IIIIDI below. 



D. Summary 



We have seen in this section that there are three important classes of behaviour for (f> outside an isolated body: 
the linear, pseudo-linear and non-linear regimes. In fact, although the mathematical analysis differs, <fi behaves in 
same way in both the pseudo-linear and linear regimes. We have shown that linear, or pseudo-linear, behaviour will 
occur whenever conditions (fTTk-c) on m c R hold. As m c R is increased, conditions (fTTb -c) will eventually fail. As m c R 
is increased still more, a thin-shell forms and we move into the non-linear regime. A thin-shell will exist whenever 
the thin-shell conditions, cqs. p^b -c). hold; these are equivalent to m c R > (m c R) e ff. We have seen that, in the 
non-linear regime, the far field is independent of the coupling of the chameleon to the isolated body. The main results 
of this section are summarised below. 

We have been concerned with a spherical body of uniform density p c and radius R. The background has density 
Pb <C p c - The chameleon in background (r ^ R) takes the value 4>b and its mass there is mj, = m^(0b). We also have 



Linear and Pseudo-Linear Behaviour 



Non-linear effects are negligible when: 

m c R « min ( (18) 1 / 6 f 



n+4 
3(n + 2) 



m c R < min V3 ^ 
m c R < 1, n = -4, 
and outside the body, r > R, (j> behaves like: 

, {m c R) 2 dp c Re- m * r 

<P ~ <Pb T7 — T 

3(n + l)r 



y/2\n+l 

( Bh,) 



\ra a ) 



, n < -4, 



V2|nTTj 



, n > 0, 



1 



2(j/ + ln(min(r/fl,l/m 6 ))) 2r 
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where yo is given by: 

(m c R) 3 = /A + I(^AV /2 
9 \ 2y 3 V 2yo) 

When non-linear effects are negligible a body will certainly not have a thin-shell. 



Bodies with thin-shells 



A body with have a thin-shell when: 



,R > 



1 



1 



m c R > m in(V3(^)™,V%TlT|(^)^ 
m c R > 4, n = -4. 

Outside the body there are two regimes of behaviour. If (r — R)/R <C 1 and m^((f>) / 'm b ^> 1 then <\> is given by 



, n < -4, 
n > 0, 



1 



\n+2\(r-R) fn+l 



to</.(0) ^2n(n + 1) 



n/2+l 



1 



and: 



If, however, (r — R)/R > 1 then 



= sgn(n)Af 



<f>~ (f> b 



n(n + l)\M 
ml{4>) 



P eff Me- m " r 



T2\ "+ 2 



47rM p/ r 

where M. is the mass of the body, and (3 e f / is the effective coupling 



Peff 
0eff 

Peff{r) 



M 



n < -4, 



^M p , ru fn(n + 1)M 2 \ "+ 2 

-JT MRl 

2nM pl 



n>0 



M 



(2Aln(min(r/i?, l/m b R))) 



-1/2 



We note that j3 e ff is independent of coupling of the chameleon to the body, and that (j> is independent of the body's 
mass, M.. When (r — R)/R> 1, <j> only depends on r, R, M , A, rt and m b . 



IV. EFFECTIVE MACROSCOPIC THEORY 



Eq. ([9]) defines the microscopic, or particle-level, field theory for 0, whereas in most cases, which we wish to study, 
we are interested in the large scale or coarse grained behaviour of (f>. In macroscopic bodies the density is actually 
strongly peaked near the nuclei of the individual atoms from which it is formed and these atoms are separated from 
each other by distances much greater than their radii. Rather than explicitly considering the microscopic structure of 
a body, it is standard practice to define an 'averaged' field theory that is valid over scales comparable to the body's 
size. If our field theory were linear, then the averaged equations would be the same as the microscopic ones e.g. as in 
Newtonian gravity. It is important to note, though, that this is very much a property of linear theories and is not in 
general true of non-linear ones. 
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FIG. 3: Illustration of the model for the microscopic structure of the body considered in this section. The constituent particles 
are assumed to be spherical of radius R, mass m p and of uniform density. They are separated by an average distance 2D. The 
average density of the body is defined to be p c . 

Non- linear effects must therefore be taken into account. Using similar methods to those that were used in section 
IIIII above, we derive an effective theory that describes the behaviour of the course-grained or macroscopic value of <j) in 
a body with thin-shell. We will identify the conditions that are required for linear theory averaging to give accurate 
results and consider what happens when non-linear effects are non-negligible. 

In this section, we derive an effective macroscopic theory appropriate for use within bodies that possess a thin-shell 
and which are made up of small particles, radius R and mass m p . These particles are separated by an average distance 
2D 3> R. The average density of the body is p c . We illustrate this set-up in FIG. [3] A thin-shell, in this sense, 
means that the average value of <fi inside a sphere of radius > D, will be approximately constant, <p ~ <p Cl everywhere 
inside the body apart from in a thin-shell close to the surface of the body. Generally the emergence of a thin-shell is 
related to a breakdown of linear theory on some level. The conditions for a body to have a thin-shell are given by 
eqs. (|19brc). The outcome of this section will be to slightly modify these conditions. Precisely, we will find that there 
is maximal, or critical, value for the average chameleon mass, m c . Oddly, this critical, macroscopic chameleon mass 
depends only on the microscopic properties of the body. 



A. Averaging in Linear Theories 

We are concerned with finding an effective theory that will give correct value of <f> c . We have defined m c = m c f > {<j) c ). 
The microscopic field equations for cj>, as given by eq. ©, are: 

-□ =-„AM=m" +, + ^> 



/ M pl 

where the microscopic matter density, p(x), is strongly peaked about the constituent particles of the macroscopic 
body but negligible in the large spaces between them. Before considering what occurs in a non-linear theory, such as 
the chameleon theories being studied here, we will review what would occur in the linear case. For the field equation 
to be linear, the potential must be, at worst, a quadratic in the scalar field <p. With the potentials considered in this 
paper, a linear theory emerges if n = 0, —1 or —2. To examine why averaging, or coarse graining, is not an issue if 
field equations are linear, and make reference to what actually occurs in chameleon theories, we shall linearise eq. 
© about <p = 4>o for some (f>o- It is important to note that we are performing this linearisation only for the purpose 
of showing what occurs in linear theories; we are not claiming that a linearisation, such as this, is actually valid. 
Defining (f> = (f>o + <j>\ , and neglecting non- linear terms, we obtain: 

- = —n\M 3 ( MX** + n(n + 1)AM 2 (M.Y** <h + ^- (25) 
\<PoJ \4>oJ M p i 
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We will write the averaged, or coarse-grained, value of a quantity Q(x) as < Q > (x) and define it by: 

<Q>m .l«gw*!=£, (26) 

J d A yb)(x - y) 

where the function Q(x — y) defines the coarse graining and the integral is over all space. Different choices of will 
result in different coarse-grainings. If we are interested in averaging over a radius of about D around the point x, 
then a sensible choice for would something like: 

£-y| 2 

0i (a; — y) — e ~~d 1 , 



or 



2 (x-g)=H(D-\x'-y\), 

where H (x) is the Heaviside function. For the coarse graining process to be well-defined we must require that, whatever 
choice one makes for 0, it vanishes sufficiently quickly as \x — y\ — > oo that the integrals of eq. (|26p converge. This 
will usually require ~ 0{\x — y\~ 3 ) as \x — y\ — > oo. 

Consider the application of the averaging procedure to the linear field equation given by eq. (|25[) . It follows from 
the assumed properties of that < Oipi >= □ < <j>i > and so 

□ < >= —nXM 3 (M_) n+1 +n (n + l)XMi ("Y" < > + f><P>& . 

V<Po/ \(paj M p i 

This is the averaged field equation for fa. It is clear from the above expression, that, although the precise definition 
of the averaging operator depends on a choice of the function 0, the averaged field equations are independent of 
this choice. This independence is a property of linear theories but it is not, in general, seen in non-linear ones. The 
averaged field equations for a non-linear theory will, generally, depend on ones choice of averaging. In this section, we 
take our averaging function to be 02 as defined above; this is equivalent to averaging by volume in a spherical region 
of radius D. It is also clear that, for a linear theory, the averaged field equation for <f>x is functionally the same as the 
microscopic equation. This, again, would not be true if non-linear terms where present in the equations; in general 

< 0™ >^< (jy >" 

unless n = or 1 or (p is a constant. 



B. Averaging in Chameleon Theories 

Our aim, in this section, is to calculate the correct value of < (f> > and < m^((f>) > inside a body with a thin-shell. 
We have defined <j> c —< 4> > and p c =< p >. Although these calculations will implicitly depend on our choice of 
averaging function, our results should also be approximately equal, at least to an order of magnitude, to those that 
would be found using any other sensible choice of coarse-graining defined over length scales of about D or greater. 

If our chameleon theories were linear, we have seen that we would expect tfi = (\>l %n ^ where 

-U4>f n) = -nXM 3 

and so, for <pi % n ^ ps const, we have: 




= M 



i\M ril M 3 



-l/(n+l) 



In appendix [Cl we show that, for some values of R, D and m p , linearised theory will give the correct value of <j) c to 

a high accuracy i.e. <j) c w (pi • This happens when there either exists a consistent, everywhere valid, linearisation 
of the field equations or we can construct a pseudo-linear approximation along the same lines as was done in section 
IIIIBI . However, for some values of R, D and m p , we find that non-linear effects are unavoidable. When R, D and 
m p take such values, we will say that we are in the non-linear regime. We find that, just as it did in section IIII CI 
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above, the non-linear regime features /3-independent critical behaviour. The details of these calculations can be found 
in appendix [U] 
We define: 



(n(n + l))^3 M" 1 



n{n + 1) 
MR 



M 



{ 3f3m p 
V47rM p/ |n| 

_ x / 3/3m p 



n + 2 



(27) 



\47rM p/ |n 



i±2 
3 



and note that D*/D c = (D c / i?)( Tl + 1 )< /3 . For the linear approximation to be valid we need both m c D <C 1 and 



\D 3 /2R < 1. This is equivalent to: 



D c < D < D„, ?i < -4 
max(D c ,Zr) <C D, n> 0. 



When n = —4 we require D <C U* and: 



ri2 s3/2 A i/2 /" 3 / 3to p 



< 1. 



We can see that, for given m p and R, it is always possible to find a D such that the linear approximation is valid 
when n > 0. However, when n < —4, it is possible that there will exist no value of D for which the above conditions 
hold. Whenever the linear approximation holds we have: 

We can construct a pseudo-linear approximation whenever: 

m m <2\n + l\(l-( 
2(n + l) (l 



D) 



n < -4 



^ > 1 and ^ < 



24:? 



1/6 



3 



n > 



n = -4. 



(28a) 
(28b) 
(28c) 



PS) 

When the pseudo-linear approximation holds we again find: 

m c » m (4«*>) . 

As the inter-particle separation, £), is decreased we will eventually reach a point where eqs. (p^Sk -c - ) fail to hold. 
When this occurs it is because non-linear effects have become important inside the individual particles that make 
up the body. As D decreases still further these particles will eventually develop thin-shells of their own. Non-linear 
effects become important when: 



H-l 



> 3|n + l| 1- 



n < -4 



D, 



> 



m<f, 



D 

D > 



3(n+l) ^1 - 

243 \ 1/6 
2\n{D/R) J 



R\~ 
D 



3 



n > 



n = -4. 



(29a) 

(29b) 
(29c) 



These conditions define the non-linear regime. Between the pseudo-linear, and fully non-linear regimes, there is, of 
course, some intermediate region, however this has proven too difficult to analyse analytically. We therefore leave 
the detailed analysis of this intermediate behaviour to a later work. This intermediate region is, however, in some 
sense small and so we do not believe it to have any great importance with respect to experimental tests of chameleon 
theories. 



22 



When the individual particles develop thin-shells, the 0-field external to the particles will be, by the results of 
section Hill independent of 0. This ensures that the chameleon mass far from the particles is also independent of (3. 
Therefore, whenever a body falls into the non-linear regime, the average chameleon mass will take a critical value, 
m c = m ™ 1 . This is defined in a similar way to which (m c R) cr i t was in section fill CI i.e. mj;™ 4 is the maximal mass 
that the chameleon may have when r ~ O(D) such that, when the microscopic field equations are integrated, (<p c /(f>) n 
is finite for all r > R. This definition implies a relationship between m c c rlt , R and D, however it does not depend on 
either M or A. m c " t is also found to depend on n; this is because n defines precisely how quickly (<fi c /4>) n blows up. 
We derive expressions m c " 1 in appendix [Cl finding: 



crit 



y/3\n+l\ ( R 



D 



q(n)/2 



= X/D, 



D / 

n = -4, 



S(n), 



n ^ -4 



(30) 



where q(n) = mm((n + 4)/(n + 1), 1), S(n) = 1 if n > and S(n) = (7(n)/3) 1 / 2 l" +1 l if n< 0. X is given by: 

3^ 



y/2\n{D/R) 



X cosh X — sinh X. 



We plot m^ rlt D vs. \n(D/R) in figure f?] For an everyday body with density similar to water, we approximate R and 
m p respectively by the radius and mass of carbon nucleus (say) and find ln(D/i?) ss 11, and so m c " t w 1A/D when 
n = -4. When n ^ -4, mf lt D < 1 follows from R/D < 1. 



Plot of m . D vs. ln(D/R) for n = 



Plotofln(m . D) vs. ln(D/R) for rut -4 





FIG. 4: Dependence of the critical chameleon mass on D/R. The above plots show how m cr i t D depends on \n{D / R) for 
different values of n. The cases n — —4 and n 7^ —4 are qualitatively different and are therefore shown on separate plots. 
m C rit is the maximal mass the chameleon can take inside a thin-shelled body. 2D is the average separation of the particles that 
comprise that body and R is the average radius of the constituent particles. Typically we find that \n{D/R) ~ 11 for bodies 
with density p ~ 1 — 10 gem -3 ; ln(D/R) = 11 is indicated on the plots. Note that when n > 0, m cr itD is independent of n. 
Also note that in <f> 4 theory m cr i t D ~ 0(1) whereas for other value of n it is generically much smaller. 

It is interesting to note that, even though m c J lt is a macroscopic quantity, it depends entirely on the details of the 
microscopic structure of the body i.e. R and D. By combining the results of this section, we find that the average 
mass of the chameleon inside a body with thin-shell that is itself made out of particles is given by: 

• f f \n\XM p iM 3 \ ^5^) crit/ n \ 
m c = mm l M [ ^ ) » ™>c (». R, D ) J • 

When evaluating the thin-shell conditions, eqs. (fT^h -c). it is therefore this value of m c that should be used. 



V. FORCE BETWEEN TWO BODIES 



In the previous two sections, we have considered how the chameleon field, </>, behaves both inside and outside an 
isolated body. In this section we study the form that the chameleon field takes when two bodies are present, and use 
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our results to calculate the resultant ^-mediated force between those bodies. The results of this sections will prove to 
be especially useful when we come to consider the constraints on chameleon field theories coming from experimental 
tests of gravity in section I VII below. 

Chameleon field theories, by their very nature, have highly non-linear field equations. This non-linear nature is 
especially important when bodies develop thin-shells. As a result of their non- linear structure, one cannot solve the 
two (or many) body problem by simply super-imposing the fields generated by two (or many) isolated bodies, as one 
would do for a linear theory. 

When the two bodies in question have thin-shells, we shall see that the formula for the </>-force is highly dependent 
on the magnitude of their separation relative to their respective sizes. We shall firstly consider the case where the 
separation between the two bodies is small compared the radius of curvature of their surfaces, and secondly look at 
the force between two distant bodies. Finally we will consider the force between a very small body and a very large 
body. We will also look at what occurs when one or both of the bodies does not have a thin-shell. 



A. Force between two nearby bodies 



We consider the force between two bodies (hereafter body one and body two) whose surfaces are separated by a 
distance d. Both bodies are assumed to satisfy the thin-shell conditions. The two bodies are taken to be nearby in 
the sense that: d <€. i?i,i?2 where R\ and R2 are respectively the radii of curvature of the surface of body one and 
body two. Since d <C R±, R2 we can ignore the curvature of the surfaces of bodies to a first approximation. With this 
simplification we treat the bodies as being infinite, flat slabs and take body one to occupy the region x < 0, and body 
two the region x > d. We use a subscript 1 to refer to quantities that are defined for body one: e.g. the density of 
body one is pi and the chameleon mass deep inside body one is mi , and a subscript 2 for quantities relating to body 
two. Additionally a subscript or superscript s is uses to refer to quantities that are defined on the surfaces of the two 
bodies e.g. mf is the chameleon mass of the surface of body one. Subscript is used two label quantities defined at 
that point between body one and body two where d<p/ dx = 0. We assume also that the background chameleon mass, 
m,b, obeys rribd <C 1, we discuss later what occurs if this is not the case. 

We now consider the ^-mediated force on body one due to body two. With the above definitions, and <p obeys: 

d2<t) V 

d^ = 



in < x < d and 



in x < 0. Integrating these equations we find: 



fd<f) s 2 



\ dx 



2(V(<f>)-V ), (31) 



in < x < d, and in x < we have: 



2 \dx J M, 
Matching these expressions at x — we have 



ifdtY = vw _ Vi + p P M-M 



(3 P i 

If the second body where not present then Vb = and <p\ = <j>\ where: 

^ A, M ? lVl 

91-91 = —5 ■ 

PPi 



The attractive force per unit area of body one due to body two is therefore: 
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This holds for all V(<j>) not just the 4> ™ potentials considered in this work. To find Vq we integrate eqn. (|3 1 1) in the 
region < x < d and find: 

j- 2 d = ( r dx + r dx \ (33 ) 

\V 0t4 ,\ Ui W{x)y/x=\ Ji W(x)Vx~^TJ ' 1 ' 

where y\ = V" /Vq and y 2 — V 2 S /Vq and W(x = V/Vo) = Vq/Vq^. We evaluate the integrals in the above expression 
in two important limits. 

1. Limit 1: yi = V{ /Va = 1 + 8 
In this limit we assume V\ ~ V{ « Vq, this would occur if V\ < V% and d is suitably small. In this limit: 



W(x)^x~l J a W(l + z 2 ) 
We also define </>f — <f>i = <^>f e > 0. With this definition equation (|32|) becomes: 



i£ = ^ 1 — ~ 0i e - 



thus 



'2Vb<5 



and 



v v 2Vl5 ±n(X\ 

Mpi V mf 



Limii 2; y x = Vf /V > 1 



This limit occurs when either d is suitably large or if V\ ^> V 2 . We take 1/yi = 6 <C 1. We consider V oc <f> n 
potentials and so W — x 1+1 / n . In this limit we have: 



r dx (\ ^-^1 ^ o 

7i iy(x)V^T Js 



1 1 1\ / 2n 
W(x)y/x~T ~ J s ^ ' ~ " V" + 2' 2 J ~ ^n + 2 

where _B( • , • ) is the Beta function. To leading order in 5 we have: 

M pi Vi 



0Pi 



We are now in a position to evaluate Vo(d) and hence F$/A. Without loss of generality we take V\ < V2 and consider 
three limits: 



3. Large separations 



If 



mid, 3> 



2 v / 2r, 



n + 2 \n+l 
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then we have V2 > Vi ^> Vq and so: 
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2X\n\M 



f+i 



d = 2B 



25 



In this limit: 



where we have defined 



^ - V = A^M 4 ( ^A±MT \ " +2 = X ^K n M\Md)^h. 
A \ \n\Md J 



K n = 



If- 



then V2 ~ Vi w Vq and: 



mid > 



v « ^ 1/2 y 2 1/2 



The force per area is therefore: 



4. Small separations: m-z ~ mi 



n + 1 \ / f m 2 

mi 



n/(n+2) 



n \ V^ 2 m d 



n + lj Vl 1/2 + V 2 1/2 



+ 1; vi /2 + v 2 1/2 
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m 2 d <C 1 



n \ V 1 ' mid 



n+lj Vl 1/2 + V 2 1/2 



n \ y/^mid 



+ u n 1/2 +^ 2 i/2 



5. Small separations: m-z S> mi 
When m 2 3> m x it is generally the case that V 2 > V{ 3> V > V\. In this case we have: 

This limit is valid provided that Vo > and V{ 3> Vo, which requires: 

1+2 



mi / n 



(n + 2)mid 



< 1. 



m 2 V n + 1 / ^/2n(n + 1) 
If m^d is not <C 1 then the F^/A given above are further suppressed by a factor of exp(— rribd). 



B. Force between two distant bodies 

We shall now consider the force between two bodies, with thin-shells, that are separated by a 'great distance', d. 
By 'great distance' we mean d 3> Ri, R2, where R\ and i? 2 are respectively the length scales of body one and body 
two. Given that d^> R\, i? 2 , then, to a good approximation, we can consider just the monopole moment of the field 
emanating from the two bodies, and model each body as a sphere with respective radii R\ and i? 2 . 

We expect that, outside some thin region close to the surface of either body, the pseudo-linear approximation (with 
the field taking its critical value i.e. (m c R) — > (m c R) cr it ~ (m c i?) e //) is appropriate to describe the field of cither 
body. In the region where pseudodinear behaviour is seen, we can safely super-impose the two one body solutions to 
find the full two body solution. 

Close to the surface of body one the mass of the chameleon induced by body one will act to attenuate the perturba- 
tion to the chameleon field created by body two. This effect can be quite difficult to model correctly. We can predict 
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the magnitude of the field, however, by noting that the perturbation to (f>, induced by body two, near to body one 
will be 

002 = — • 

From the results of section UlI C[ we know that fl^ffM-i/Mpi only depends on the radius of body two and the theory 
dependent parameters, M, A and n; e ff is as given in eqs. (|24h-c); raj, is the mass of the chameleon field in the 
background i.e. far away from cither of the two bodies. M.2 is the mass of body two. 

We define the perturbation to (f> induced by body one near to body one in a similar manner 

P^ f M ie - m » d 



AirMpid 

From eqs. (|24a -c) we know that 84>\ is independent of j3 and the mass of body one, M.j_. The force on body one due 
to body two will be proportional to V<5</>2, however, since this must also be the force on body two due to body one, it 
must also be proportional to V5</>i evaluated near body two. ^From this we can see that the force on one body due 
to the other must, up to a possible 0(1) factor, be given by 

_ ^^1^^2(1 + m b d)e~ m > d 
(47rM pl ) 2 d 2 ■ 

The functional dependence of this force on M, n, Ri, R2 and A depends on whether n < —4, n = —4 or n > 0. We 
consider these three cases separately below. In all cases the force is found to be independent of j3, Mi and M.2- 



When n < -4, eq. ([24k) gives 



Case n < —4 



ith 1 - 

and i?2, separated by a distance d ^> Ri, i? 2 , is therefore given by 

V 
d 2 
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The expression for f3 e ^ is similar but with 1 — > 2. The force between two spherical bodies, with respective radii R\ 

7 (n)\^ {MRx)^{MR 2 )^{l + m b d)e-^ d 
t 4> = -JZT m • i 34 ) 



2. Case n > 

eff " "' '"' "" 



If n > then fity. is given by eq. l[2"lb) to be 



AnM ri 



= MRx 



(2) 

Again the expression for f3 e ^ is similar. The force between two distant bodies is therefore found to be 

fn(n+ 1)M 2 \^ M 2 R 1 R 2 (l + m b d)er m » d 
mi J TP ■ (35) 



3. Case n = — 4 

?(1,2) 



When n = —4, P e ff are given by eq. ([Mb ) and they are actually weakly dependent on r. Using eq. ([2~4"b) we find 
that 

= (l + m b d)e-^ d 
8A A /ln(d/i?i)ln(d/i? 2 )ef 2 ' 
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C. Force between a large body and a small body 

One subcase that is not included in the above results is the force between a very large body with radius of 
curvature Ri, and a very small body with radius of curvature R2, that are separated by an intermediate distance d 
i.e. R2 ^> d ^> Ri. We assume that both bodies have thin-shells. In this case we find a behaviour that is half-way 
between the two cases described above in sections IV Al and IV Bl The magnitude of field produced by the large body 
will be much greater than that of the small body. 

If we ignore the small body and assume that the average mass of the chameleon in the background obeys m b d <C 1, 
then the field produced by body one is given by eq. (f2"2")l . Using this equation, we find that 



dd> , , 1 / 2 \ 2 <™+ 2 > , "+4 o 

The effective coupling of body two to this ^-gradient will be j3 e ft as it is given by eqs. (f2"lk -c). If m b d > 1 then this 
gradient will be, up to an order 0(1) coefficient, attenuated by a factor of (1 + m b d) exp(— m b d). The force between 
the two bodies is therefore given by 

F <t> = (Mi? 2 )"+ 2 ' 
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As before, d is the distance of separation. These formulae will prove useful when we consider the 0-force between 
the Earth and a test-mass in laboratory tests for WEP violation in section IVI CI 

D. Force between bodies without thin-shells 

If neither of the two bodies have thin-shells, 4> behaves just like a standard, linear, scalar field with mass m b . The 
force between the two bodies, with masses M.\ and M2, is given by: 

_ p 2 MiM2(l + m b d)e- m " d 

As above, m b is the mass of the chameleon in the background. If one of the bodies has a thin-shell, body one say, but 
the other body does not, then the force is given by 

$f f PMiM*(\ + m h d)e- m * d 

Fa, — 



whenever d^> R\ where R\ is the radius of curvature of body one and (f^)* is as given by eqs. (pM)) . If <i <gc i?i then 
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Fa, = A"T5 ( 2 ^ 2< " +2) (Md)^ ^2(1 + mb d)e-™>° 



As above, d is the distance of separation. 

E. Summary 



In this section we have considered the force that the chameleon field, 0, induces between two bodies, with masses 
Ai\ and M.2 and radii R\ and R2, separated by a distance d. The chameleon mass in the far background is taken to 
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be m b . When both bodies have thin-shells, we found that, to leading order, the 0-force between them is independent 
of the matter coupling, f3, provided that mid, m 2 d 3> 1; mi is the mass that the chameleon has inside body one, and 
m 2 is similarly defined with respect to body two. The force between two such bodies is also independent Mi and M 2 
but does in general depend on i?i and R 2l as well as on M, n and A. The main results of this section are summarised 
below. 



Neither body has a thin-shell 



P 2 MiM 2 {l + m b d)e- m » d 
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Body one has a thin-shell, body two does not 
If the two bodies are close together (m^ 1 <C d <C Ri) then 
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If the bodies are far apart d^> Ri then 
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This force is independent of the mass of body one, and the chameleon's coupling to it. 



Both bodies have thin-shells 



If the two bodies of close by and m-^ 1 ,m 2 1 <C d -C Ri, R2 then the force per unit area is 
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Different formulae for Fa,/ A apply if either m\d < 1 or m 2 d < 1, these are given at the end of section fV Al above. 
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where 7(71) is given in tabled] If body one is much larger than body two and they are at an intermediate separation: 
R 2 , m' 1 > d > i?i then 

F " - Vh) {MR2) Ifw 

[Md) n + 2 -5 , n < -4 

er 
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VI. LABORATORY CONSTRAINTS 



The best bounds on corrections to General Relativity come from laboratory experiments such as the Eot-Wash 
experiment, [38| and Lunar Laser Ranging tests for WEP violations, [26|, [27| ■ At very small distances d < 10/im, the 
best bounds on the strength of any fifth force come from measurements of the Casimir force. 

In this section we will consider, to what extent, the results of these tests constrain the class of chameleon field 
theories considered here. We will find the rather startling result that (3 3> 1 is not ruled out for chameleon theories. 

One of original reasons for studying chameleon theories with n > potentials, [25|, was that the condition for an 
object to have a thin-shell, eq. (|19b ). was found to depend on the background density of matter. It is clear from 
eq. (|19b ). that the smaller pb is, the larger m c R must be for a body to have a thin-shell. This property lead the 
authors of ref. [25[ to conclude that the thin-shell suppression of the fifth-force associated with cf> would be weaker 
for tests performed in the low density vacuum of space pb ~ 10 _25 gcm -3 , than it is in the relatively high density 
laboratory vacuum pb ~ 10~ 17 gcm~ 3 . As a result, it is pos sible that, if the same experimental searches for WEP 
violation, which were performed in the laboratory in [HoL l5ll. l62l l63l |64| . were to be repeated in space, they would 
find equivalence principle violation at a level greater than that already ruled out by the laboratory-based tests. It is 
important to note that this is very much a property of n > theories. It is clear from eqs. (TH^i) and (fTTjb ) that, when 
n < —4, the thin-shell condition, for a body of density p c , is only very weakly dependent on the background density of 
matter when pb <C p c - As a result, space-based searches for WEP violation will not detect any violation at a level that 
is already ruled by lab-based tests for n < —4 theories. Planned space-based tests such as STEP [54], SEE [53], GG 
[55I ] and MICROSCOPE [5(| promise a much greater precisions than their lab-based counterparts. MICROSCOPE 
is due to be launched in 2007. This improved precision will, in all cases, provide us with better bounds on chameleon 
theories. 



A. Eot-Wash experiment 

The University of Washington's Eot-Wash experiment, [H, [3§], is designed to search for deviations from the 1/r 2 
drop-off gravity predicted by General Relativity. The experiment uses a rotating torsion balance to measure the 
torque on a pendulum. The torque on the pendulum is induced by an attractor which rotates with a frequency u>. 
The attractor has 42 equally spaced holes, or 'missing masses', bored into it. As a result, any torque on the pendulum, 
which is produced by the attractor, will have a characteristic frequency which is some integer multiple of 21lu. This 
characteristic frequency allows any torque due to background forces to be identified in a straightforward manner. The 
torsion balance is configured so as to factor out any background forces. The attractor is manufactured so that, if 
gravity drops off as 1/r 2 , the torque on the pendulum vanishes. 

The experiment has been run with different separations between the pendulum and attractor. The Eot-Wash 
group recently announced some new results which go a long way towards better constraining the parameter space 
of chameleon theory, [39( | . The experiment has been run for separations, 55/im < d < 9.53mm. Both the attractor 
and the pendulum are made out of molybdenum with a density of about pub ~ 10 gem™ 3 and are 0.997mm thick. 
Electrostatic forces are shielded by placing a 10 /im thick, uniform BeCu sheet between the attractor and pendulum. 
The density of this sheet is pBeCu ~ 8.4 gem™ 3 . The role played by this sheet is crucial when testing for chameleon 
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fields. If (3 is large enough, the sheet will itself develop a thin-shell. When this occurs the effect of the sheet is not 
only to shield electrostatic forces, but also to block any chameleon force originating from the attractor. 

The force per unit area between the attractor and pendulum plates due to a scalar field with matter coupling (3 
and constant mass to, where 1/m <C 0.997toto is: 

F ae- md 2nGp 2 Mb 
l A l to 2 

where a — 1 /Air and d is the separation of the two plates. The strongest bound on a coming from the Eot-Wash 

experiment is a < 2.5 x 10~ 3 for 1/to = 0.4 0.8toto. 

Depending on the values of f3, M and A there are three possible situations: 

• The pendulum and the attractor have thin-shells, but the BeCu sheet does not 

• The pendulum, the attractor and the BeCu sheet all have thin-shells. 

• Neither the test masses nor the BeCu sheet have thin-shells. 

In the first case the ^-mediated force per unit area in a perfect vacuum is given by one of the equations derived in 
section IV Al (depending on the separation d) . In reality the vacuum used in these experiments is not perfect actually 
has a pressure of 10~ 6 Torr which means that the chameleon mass in the background, nib, is non-zero and so F^/A is 
suppressed by a factor of exp(— m^d). Fortunately however mbd <C 1 for all but the largest /3. A further, and far more 
important, suppression occurs when the BeCu sheet has a thin-shell. If mseCu is the chameleon mass inside the BeCu 
sheet and dseCu the sheet's thickness, the existence of a thin-shell in the electromagnetic shield causes the chameleon 
mediated force between the pendulum and attractor to be suppressed by a further factor of exp(— niBeCudBeCu)- The 
thin-shell condition for the BeCu sheet implies that niBeCudBeCu 3> 1, and so this suppression all but removes any 
detectable chameleon induced torque on the pendulum due to the attractor. The BeCu sheet will itself produce a 
force on pendulum but, since the sheet is uniform, this force will result in no detectable torque. If M and A take 
natural values the electromagnetic sheet develops a thin-shell for (3 > 10 4 : as a result of this the Eot-Wash experiment 
can only very weakly constrain large /3 theories. 

If neither the pendulum or the attractor have thin-shells then we must have m^d <C 1 and the chameleon force is 
just (3 2 /Air times the gravitational one. Since this force drops of as 1/r 2 it will be undetectable from the point of 
view this experiment. In this case, however, (3 is constrained by other experiments such as those that search for WEP 
violation (2g, H3, [62|, [63|, HH or those that look for Yukawa forces with larger ranges [65J . 

Given all of the considerations mentioned above, we used the formulae given in section IV Al to evaluate the latest 
Eot-Wash constraints on the parameter space of chameleon theories. Our results are shown in FIG In these plots 
the shaded region is allowed by the current bounds. 

When (3 is small, the chameleon mechanism present in these theories becomes very weak, and from the point of 
view of the Eot-Wash experiment, <fi behaves like a normal (non-chameleon) scalar field. When (3 ^> 1, the </>-force 
is independent of the coupling of the chameleon to attractor or the pendulum, but does depend on the mass of the 
chameleon in the BeCu sheet, rriBeCu- The larger mseCu is, the weaker the Eot-Wash constraint becomes. Larger (3 
implies a larger mseCm an d this is why the allowed region of parameter space increases as (3 grows to be very large. 

When n ~ —4, we can see that a natural value of A is ruled out for 10 _1 < (3 < 10 4 , but is permissible for (3 > 10 4 . 
This is entirely due to the that BeCu sheet has a thin-shell, in n = —4 theories with A = 1/4!, whenever (3 > 10 4 . 

It is important to stress that, despite the fact that the Eot-Wash experiment is currently unable to detect (3 3> 1, 
this is not due to a lack of precision. One pleasant feature, of the /3-independence of the </>-force, is that if you can 
detect, or rule out, such a force for one value of f3 3> 1, then you will be able to detect it, or rule out, all such 
(3 3> 1 theories. If design of the experiment can be altered so that electrostatic forces are compensated without using 
a thin-sheet then the experimental precision already exists to detect, or rule out, almost all (3 3> 1, <j) A theories with 
Awl/4!. 

In conclusion, an experiment, along the same lines of the Eot-Wash test, could detect, or rule out, the existence of 
sub-Planckian, chameleon fields with natural values of M and A in the near future, provided it is designed to do so. 

B. Casimir force experiments 

Short distance tests of gravity fail to constrain strongly coupled chameleon theories as a result of their use of a 
thin metallic sheet to shield electrostatic forces. However, experiments designed to detect the Casimir force between 
two objects, control electrostatic effects by inducing an electrostatic potential difference between the two test bodies. 
By varying this potential difference and measuring the force between the test masses, it is possible to factor out 
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FIG. 5: [Colour Online] Constraints on chameleon theories coming from the Eot-Wash bounds on deviations from Newton's 
law. The shaded area shows the regions of parameter space that are allowed by the current data. The solid black lines indicate 
the cases where M and A take 'natural values'. The dotted-black line indicates when M = M§ := M p i//3 i.e. when the mass 
scale of the potential is the same as that of the matter coupling. Other n < —4 theories are similar to the n = —8 case, whilst 
the n — 4 plot is typical of what is allowed for n > theories. The amount of allowed parameter space increases with \n\. 

electrostatic effects. As a result, Casimir force experiments provide an excellent way in which to bound chameleon 
fields where the scalar field is strongly coupled to matter. 

Casimir force experiments measure the force per unit area between two test masses separated by a distance d. It is 
generally the case that d is small compared to the curvature of the surface of the two bodies and so the test masses 
can be modeled, to a good approximation, as flat plates. 

In section IV Al we evaluated the force per unit area between two fiat, thin-shelled slabs with densities p\ and p2- 
The Casimir force is between two such plates is: 

Fca. = 7T 2 

A 240d 4 ' 

Whilst a number of experimental measurement of the Casimir force between two plates have been made, the most 
accurate measurements of the Casimir force have been made using one sphere and one slab as the test bodies. The 
sphere is manufactured so that its radius of curvature, R, is much larger than the minimal distance of separation d. 
In this case the total Casimir force between the test masses is: 

(\ ix 2 1 \ (R (amf\ , 

In all cases, apart from when n = —4 and mid, m-id S> 1, the chameleon force per area grows more slowly than d~ 4 as 
d — » 0. When n — —4 and mid, m<zd > 1 we have F^/A cx <i~ 4 . It follows that the larger the separation, d, used, the 
better Casimir force searches constrain chameleon theories. Additionally, these tests provide the best bounds when 
the test masses do have thin-shells as this results in a strongly d dependent chameleon force. Large test masses are 
therefore preferable to small ones. 
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Note that if the background chameleon mass is large enough that rribd > 1 then is suppressed by a factor of 
exp(— mjjd). The smaller the background density, pi,, is, the smaller become. Since small nib is clearly preferably, 
the best bounds come from experiments that use the lowest pressure laboratory vacuum. 



Casimir force bounds on n = -8 theories 



Casimir force bounds on n = -4 theories 





Casimir force bounds on n = 4 theories 




FIG. 6: [Colour Online] Constraints on chameleon theories coming from experimental searches for the Casimir force. The 
shaded area shows the regions of parameter space that are allowed by the current data. The solid black lines indicate the cases 
where M and A take 'natural values'. The dotted-black line indicates when M = M,/, := M p i//3 i.e. when the mass scale of 
the potential is the same as that of the matter coupling. Other n < —4 theories are similar to the n = —8 case, whilst the 
n — 4 plot is typical of what is allowed for n > theories. The amount of allowed parameter space increases with |n| i.e. as 
the potential becomes steeper. 

In (40| . Lamoreaux reported the measurement of the Casimir force using a torsion balance between a sphere, with 
radius of curvature 12.5cm ± 0.3cm and diameter of 4cm, and a flat plate. The plate was 0.5cm thick and 2.54cm in 
diameter. The apparatus of placed in a vacuum with a pressure of 10 _4 mbar. Distances of separation of 6 — 60/im 
we used and in the region d w 7 — 10/im it was found that: 

^measured _ ptheory | < 

Another measurement, this time using a microelectromechanical torsional oscillator, was performed by Decca et 
al. and reported in [4l| . In this experiment, the sphere was much smaller than that used by Lamoreaux, being 
on 296 ± 2pm in radius; the plate was made of 3.5/im thick, 500 x 500^m 2 polysilicon. The smallness of these 
test masses means that they will only have thin-shells when (3 is very large. In the region d as 400ran — 1200nm, 
^pmeasured _ Fq^ v \ < 7.5 x 10 _2 /idyn. We show how these experiments constrain the parameter space of <\T n 
chameleon theories in FIG E3 Other Casimir force tests (e.g. [42[ ) are less suited to constraining chameleon theories 
such as those considered here. As in the previous plots, the shaded area is allowed, the solid black line is M ~ (pa) 1 ^ 4 
or A = 1/4!, and the dotted black line is M = M p i//3. We note that Casimir force experiments provide very tight 
bounds on A and M when ,8 > 1. A natural value of A when n — —4 is ruled out for all (3 > 10 4 . When this is 
combined with the latest Eot-Wash data we can rule out A = 1/4! for all /3 > 1. For other n, we see that we cannot 
have M much larger than its 'natural' value (pa) 1//4 for large j3. If the bounds on extra forces at d ~ 1 — IQpm can 
be tighten by roughly an order of magnitude then a natural value of M can be ruled out for all large j3. Casimir force 
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tests provided by far the best bounds on M and A for large /3. We note that n = —4 theories are the most tightly 
constrained by Casimir force experiments, this is not surprising since F^/Fc as ~ const and is > 1 in the region when 
m\d,mid 3> 1 in this model, whereas in all the other theories F^/Fcas decreases as d is made smaller. 

More generally, the steeper the potential in a given theory is, the more slowly F^/A increases as d — > and, as a 
result, the weaker Casimir force bounds on the theories parameter space are. 



C. WEP violation experiments 

The weak equivalence principle (WEP) is the statement that the (effective) gravitational and inertial mass of 
a body are equal. If it is violated then the either the strength of gravitational force on a body depends on its 
composition, or there is a composition dependent 'fifth-force'. Since we believe that gravity is geometric in nature, 
most commentators, ourselves included, would tend to interpret any detection of a violation of WEP in terms of 
the latter option. The existence of light scalar fields that couple to matter usually results in WEP violations. As 
we mentioned in the introduction to this thesis, the experimental bounds on WEP violation are exceeding strong, 
f2ot l27l l50l l5ll l52l . lo2l . lo3l . [Hj], and, at present, they represent th e str ongest bounds on the parameters of, non- 
chameleon, scalar-tensor theories such as Brans-Dicke theory, [ID, Il05l |. A number of planned satellite missions 
promise to increase the precision to which we can detect violations of WEP by between 2 to 5 orders of magnitude 
|53l. I&il. l55l. [5^| . The precision that is achievable in laboratory based tests also continues to increase at a steady rate. 

Experiments that search for violations of the weak equivalence principle generally fall into two categories: laboratory 
based experiments, which often employ a modified torsion balance, (62l. l63l l&il]. and solar system tests such as lunar 
laser ranging (LLR) [26|, [27| • 

The laboratory based searches use a modified version of the Eot-Wash experiment mentioned above. In these 
experiments the test masses are manufactured to have different compositions. The aim is then to detect, and measure, 
any difference in the acceleration of test-masses towards an attractor, which is usually the Earth, the Sun or the Moon. 
In some versions of the experiment a laboratory body is used as the attractor. 

If the test masses have thin-shells then the 0- force pulling them towards the Earth is given by equations ([37^,-c). 
The chameleon force towards the Moon or Sun is given by eqs. (j34"|) - (l36l) . If the attractor is a laboratory body then, 
depending on the separations used, the force is given by either eqs. (|34"|) -(l3l) f or by eq. (jV A[) . 

We label the mass and radius of the attractor by M.3 and R3 respectively, and take the mass and radius (or size) 
of the two test-masses to be given by {M.\ 1 Ri} and {-M 2 , ^2}- We define «i3 to be the relative strength, compared 
to gravity, of 0-force between the attractor (body three) and the first of the test masses (body one). 023 is defined 
similarly as a measure of the 0-force between body two and body three. The difference between the acceleration of 
the two test masses towards the attractor is quantified by the Eotvos parameter, ?y, where 

2|ai 3 -a 2 3| 1 1 
11 = WZ T 1 ~ I"" " a 2sl- 

M + a 13 + "23 | 

When the test-bodies have thin-shells, we found, in section [V] that the 0-force is independent of the masses of the 
test-bodies, the mass of the attractor and the coupling of the test-masses and attractor to the chameleon. The only 
property of the attractor and test bodies, which the </>-force does depend on, is their respective radii. Since the 
gravitational force between the test-masses and the attractor does depend on the masses of the bodies, it follows that 
ai3 only depends on A^i, M3, Ri, R3, M (or A), mi, and n, where m& is the chameleon mass in the background. It 
does not depend on the chameleon's coupling to the test-mass, (3. The situation with 023 is very similar. 

Since the 0-force is independent of the coupling, f3, any microscopic composition dependence in (3 will be hidden on 
macroscopic length scales. The only 'composition' dependence in a 13 is through the masses of the bodies and their 
dimensions (i?i and -R3). 

Taking the third body to be the Earth, the Sun or the Moon, experimental searches for WEP violations have, to 
date, found that n 5, 10~ 13 [H, [H, H3|- Future satellite tests promise to be able to detect violations of WEP at 
between the 10 -15 , and the 10~ 18 level, (54[. It also is claimed that future laboratory tests will be able to see 
r\ ~ 10~ 15 , however, whilst the precision to detect at such a level is achievable, there are a number of systematic 
effects that need to be compensated for, before an accurate measurement can be made. 

In most of these searches, although the composition of the test-masses is different, they are manufactured to have 
the same mass (AA\ = M.2) and the same size (R± — R2). Therefore, if the test-masses have thin-shells, we will have 
0113 = «23 and so n = identically. As a result, a chameleon field will produce no detectable WEP violation in these 
experiments. The only implicit dependence of this result on (3 is that, the larger the coupling is, the more likely it is 
that the test-masses will satisfy the thin-shells conditions. 

If one wishes to detect the chameleon using WEP violations searches, then one must either ensure that test-masses 
do not satisfy the thin-shell conditions, or that they have different masses and/or dimensions. 
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We have just argued that all of the chameleon theories considered here will automatically satisfy all laboratory 
bounds on WEP violation, provided the test-masses have thin-shells. This occurs entirely as a result of the design of 
those experiments. 

Let us consider then a putative experiment, which could be conducted, that would, in principle, be able to detect 
the chameleon through a violation of WEP. In this experiment the test masses are of different densities (p\ and p 2 ) 
but of the same mass, M.test- Crucially the radii (size) of the two bodies are taken to be different: R\ and R 2 . We 
now calculate the Eotvos parameter, 77, taking the attractor to be either the Earth, the Sun or the Moon. 



Attractor is the Earth 



If the attractor is the Earth then we obtain 
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where m is the mass of the chameleon in the background region between the test masses and the surface of the Earth; 
d is the distance between the test masses and the surface of the Earth. Ai Earth is the mass of the Earth, and REarth 
its radius. 

Current experimental precision bounds 77 < 10~ 13 . We shall assume that our putative experiment, if conducted, 
would find r] < 10~ 13 . However, even if this is the case, we are still only able to recover very weak bounds on 
{P,M, A}. The bounds on j3 are especially weak due to the /3-independence of the 0- force whenever the test-bodies 
have thin-shells. The only real bound on (3 comes out of requirement that it be large enough for the test-masses to 
have thin-shells. 

For definitcness, we take the test-masses to be spherical, with a mass of Attest = 10 g. We assume that one of them 
is made entirely of copper and the other from aluminum. If the test-masses have thin-shells then, even if m d <C 1, 
finding 77 < 1CP 13 would, when n = —4, only limit 



A > 10 
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For theories with n < —4, 77 < 10 is easily satisfied provided that 



When n > 0, the resultant 



bound is on a combination of M and m^. The WEP bounds on the parameter space of n > theories are generally 
stronger than those for other n. We plot the effect of these bounds on the parameter space of our chameleon theories 
in FIG [3 



2. Attractor is the Sun or the Moon 



Constraints on chameleon theories can also be found by considering the differential acceleration of the test masses 
towards the Moon or the Sun, rather than towards the Earth. The analysis for both of these scenarios proceeds along 
the much same lines. Since, for the reasons we explain below, the lunar bound will be by far the stronger, we will 
only explicitly consider the case where the third body is the Moon. In this case, the force between the test mass and 
the Moon is given by equations (f3~4"]) . ([33)1 and (|33|) for 71 < —4, n > and n = —4 respectively. 

We define mj, to be the average chameleon mass in the region between the Earth and the Moon, and m a tm to be the 
mass of the chameleon in the Earth's atmosphere; m; a h is the background mass of the chameleon in the laboratory. 
R a is the thickness of the Earth's atmosphere, d is the distance of separation between the laboratory apparatus and 
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the Moon. Evaluating r\ we find: 
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(38) 



(39) 



n > 0, 
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where Mmoou is the mass of the Moon and Rmoou its radius. If the Sun is used as the attractor then rj is given by a 
similar set of equations to those shown above only with M.Moon — ► -Msun and Ruoon — ► Rsun- Msun and Rsun are 
respectively the mass and the radius of the Sun. If the attractor is taken to have fixed density, then we can see that 
r\ dies off at least as quickly as 1/R 2 , where R is the radius of the attractor. It follows that the predicted value rj, 
induced by chameleon, is always smaller when the Sun is used as the attractor than when the Moon is. This is because 
Rsun ^ Rmoou- The predicted values of rj have a similar dependence on the radii of the test-masses i.e. dying off at 
least as quickly as 1/R 2 . The corollary of this result is that if we are unable to detect </> in lab-based, micro-gravity 
experiments where the radii of the test-masses and the attractor are both of the order of 10 cm, or smaller, then the 
</>-force between larger (say human-sized) objects would also be undetectably small. For this reason, in the context of 
chameleon theories, measurements of the differential acceleration of the Earth and Moon towards the Sun, e.g. lunar 
laser ranging [H, H3j) are not competitive with the bounds on WEP violation found in laboratory-based searches. 




3. Summary 

We show how bounding the Eotvos parameter by rj < 10~ 13 , with the attracting body being either the Moon or the 
Earth, constrains chameleon theories in FIG. [7] As with the Eot-Wash plots: the whole of the shaded area is currently 
allowed, whilst the more lightly shaded area is that which could be detected by proposed space-based tests of gravity 
such as SEE, STEP, GG and MICROSCOPE. It is claimed, [H, [H [H, [H[ , that these experiments will be able to 
detect rj down to 10 -18 . This improved precision is responsible for most of the increased ability of spaced based tests 
to detect a chameleon field. When n > 0, the thin-sell condition is stronger in the low density background of space 
than it is in the relatively higher density background of the laboratory; this effect accounts for some of extra ability 
that future space-based experiments have to detect the chameleon. We note that WEP violation searches only have 
any real hope of detecting the chameleon field, if M take a 'natural' value i.e. M ~ (0.1mm) -1 , in n > theories. 



D. Discussion 



In FIG. [7) we plot how all of the bounds on WEP violation mentioned in refs. [26|, H3, H3, E3, HH, as well as the 
putative bound resulting from the modified WEP violation test we considered above, constrain the parameter space 
of our chameleon theories. The Eot-Wash bounds are shown in FIG. [5j and the Casimir bounds in FIG. [6] These 
Eot-Wash, Casimir and WEP violation bounds are also included in the plots of section Hxl 

It is important to note that, the larger (3 is, the stronger the chameleon mechanism becomes. A strong chameleon 
mechanism results in larger chameleon masses, and larger chameleon masses in turn result in weaker chameleon- 
mediated forces. A stronger chameleon mechanism also increases the likelihood of the test masses, used in these 
experiments, having thin-shells. Large values of f3 cannot therefore be detected at present by the Eot-Wash and WEP 
tests, if A and M take natural values. If the matter coupling is very small, /3< 1, then the chameleon mechanism is 
very weak and 4> behaves as a normal (non-chameleon) scalar field. The current precision of laboratory tests of gravity 
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FIG. 7: [Colour Online] Constraints on chameleon theories coming from WEP violation searches. The whole of shaded area 
shows the regions of parameter space that are allowed by the current data. Future space-based tests could detect the more 
lightly shaded region. The solid black lines indicate the cases where M and A take 'natural values'. The dotted-black line 
indicates when M = := M p i/f3 i.e. when the mass scale of the potential is the same as that of the matter coupling. Other 
n < —4 theories are similar to the n — —8 case, whilst the n = 4 plot is typical of what is allowed for n > theories. The 
amount of allowed parameter space increases with \n\. 



prevent them from seeing even non-chameleon theories with (3 < 1CP 5 . Experiments that search for the Casimir force 
are better able to detect large (3 theories primarily due to the way in which they cancel electrostatic forces. The 
(3 ~ O(l) region is, however, currently inaccessible to Casimir force experiments. 

Casimir force experiments provide upper bounds on 1/A and M but do not directly rule out large [3. Upper bounds on 
j3 do arise, however, from astrophysical considerations. It so happens that for (3 < 10 20 , these astrophysical constraints 
are weaker than those coming from Casimir force searches, however they are important since they constrain how the 
chameleon field can behave in much lower (and higher) density backgrounds than those that are easily accessible in 
the laboratory, and thus effectively probe a different region of the chameleon potential. These bounds are discussed 
in sections I VIII and IVIIII below. 

In conclusion: contrary to most expectations laboratory tests of gravity do not rule out scalar field theories with 
a large matter couplings, [3 ^> 1, provided that they have a strong enough chameleon mechanism. When n = —4, a 
natural value of A is ruled out for all a = 1 /4tt < 1CP 2 , and even A ~ 0(1) requires (3 ~ 10 7 , i.e a ~ 10 13 . When 
n 7^ —4, large (3 theories with natural values of M are allowed for all n ^ —6. Chameleon theories, with natural values 
of M, could well be detected or ruled out by a number of future experiments provided they are properly designed to 
do so. 



VII. IMPLICATIONS FOR COMPACT BODIES 



In this section we will consider the effect that the fifth force associated with the chameleon field has upon the physics 
of compact bodies such as white dwarfs and neutron stars. In the preceding analysis, we have shown that the <fi- force is 
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only comparable in strength to gravity over very small scales. In neutron stars and white dwarfs, however, the average 
inter-particle separations are very small, about 10 _13 cm and 10 _10 cm respectively, and the physics of such compact 
objects can therefore be very sensitive to additional forces that only become important on small scales. The stability 
of white dwarf and neutron stars involves a delicate balancing act between the degeneracy of, respectively, electrons 
or neutrons and the effect of gravity. If the presence of the chameleon field were to alter this balance significantly, one 
might find oneself predicting that such compact objects are always unstable. If chameleon theories were to make such 
a prediction, for some values of the parameters {/3, M, A}, then we could, obviously, rule out those parameter choices. 
As well as the issue of stability, we must also consider potential astrophysical observables such as the Chandrasekar 
mass limit, and the mass-radius relationship. 

In 1930, Chandrasekar made the important discovery that white dwarfs had a maximum mass ~ 1.4M Q , [102L 1 103| ] . 
The precise valu e depends on the composition of the star. A similar maximum mass was found for neutron stars 
by Landau jl04| . It was additionally noted that the mass, M stari of a white dwarf or neutron star would depend 
on its radius, R, in a very special manner. This is the m ass-radiu s relationship. It is possible to extrapolate both 
M s tar and R from astronomical data, for example see [99l . llOdl . Il0l| . In all of those works, and others like them, the 
mass-radius relationship, as predicted by General Relativity , has been found to be in good agreement with the data. 
It is, therefore, important that the addition of a chameleon field should not greatly alter this relationship. The quoted 
lcr error bars on most of the determinations of M star and R are between about 3 and 10% of the central value. It 
would certainly be fair to say then that any new theory which predicts deviations in the value of M star (R) from the 
GR value of less than about 10%, is consistent with all current data. Much greater deviations, are however ruled out. 
We shall use this criterion to bound the parameters of our chameleon theories. 

We firstly consider how the presence of a chameleon alters the Chandrasekar mass limit and the mass-radius relation 
for both white dwarfs and neutron stars. Our analysis proceeds along the same lines as that presented in (57j . We 
will start by considering a white dwarf and then note how our results carry over to neutron stars. We suppose that 
the white dwarf contains N electrons. Charge neutrality then implies that there are N protons. There will also be 
neutrons present (approximately N of them) but for this calculation we will merely group the protons and neutrons 
together into N nucleons where each nucleon contains one proton. We denote the mass of a nucleon by m u and 
take it to be the atomic mass unit, m u = 1.661 x 10~ 24 g. White dwarfs are kept from collapsing by the pressure 
of degenerate electrons, whilst their gravitational potential comes almost entirely from the nucleons (as m u 3> m e ). 
We model the white dwarf as being at zero temperature. If should be noted that we are not interested, so much, in 
accurately determining the mass-radius relationship or Chandrasekar mass limit, as we are in seeing to what extent 
they deviate from the general relativistic form. 

In the limiting cases of non-relativistic, T = 5/3, and relativistic, T — 4/3, behaviour, the equation of state for the 
electrons can be written in polytropic form: 



with K a constant. For relativistic electrons: 



K 



P = Kp l 

31/3^2/3 



4/3 5/3' 
rriu n e ' 



where fi e ss 2 is the chemical potential for the electrons in the white-dwarf, [57( . We require that the white dwarf be 
in hydrostatic equilibrium. Ignoring general relativistic corrections, this implies 

VP = -pVcj) N — - , 

Mpi 

where the last term is the additional clement that comes from the chameleon field. 4>n is the Newtonian gravitational 
potential: V 2 ^at = Airp/M^. In most realistic scenarios the density inside the white dwarf will only change very 
slowly over length scales comparable to the inverse chameleon mass. We can therefore take 



M 



M p iM 3 n\ 



to hold inside the white dwarf. It is standard practice, [57U58||. to solve for hydrostatic equilibrium by minimizing an 
appropriately chosen energy functional. In the absence of a chameleon field this is: 



E = U + W, 
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where W is the gravitational potential energy: 



IT" = - / d 3 x~(j> N p, 



and U is the internal energy of the white dwarf. It is shown in 57] that when one perturbs the density in such a way 
that the total mass is conserved (Sp = —V ■ (p£) for some vector field £) we have: 

SE = SU + SW = J (VP + pV(j)N) ■ £d 3 x. 

It follows that SE vanishes for a star in hydrostatic equilibrium. To solve for hydrostatic equilibrium in the presence 
of a chameleon field, we need to minimise the following energy functional: 

E = E + Wa, = U + W + Wa,, 

where: 

Wa, = / d 3 x—p. 

n J Mpi 

To see that this is the correct expression we consider SWa,: 

n + l f , 3 p5{ct>p) _ f , 3 [3<f> 



SW^ = — / d 6 x^P^ = / d 6 x^-S P 
n J M p i J Mp, 

With this definition SE = is seen to be equivalent to requiring hydrostatic equilibrium: 

se = J (vp + pV4> N + • £d 3 x = 0. 

A. The Mass-Radius Relation 

Schematically we have W<f, oc 7 - L ^-{(3< <fi >/M p i)M s t ar where < .. > indicates an average and M sta r = m u N is the 
mass of the star. We note that W<j, ~ and that it is negative for n < — 4 and positive for n > 0. To study 

the effect of the chameleon upon the Chandrasekar mass limit, and the mass-radius relationship, we assume that the 
density of the white dwarf is uniform. Whilst this is not at all accurate, it is sufficient to see when the chameleon does, 
or does not, have a significant effect. Also, whilst not being accurate, this approximation still gives the mass-limit 
and mass-radius relationship up to an 0(1) factor. This said, we shall consider a more accurate model later when we 
look at general relativistic corrections. The total internal energy of the white dwarf is given by: 

U = N ((p 2 F + ml) 1/2 - Nm e ) > 0, 

where pf = N 1 / 3 /R is the Fermi momentum of degenerate electrons, [13, [Hj]. For W we find: 

3ml N 2 



W = 

5M*R 



Evaluating Wa, in this approximation yields: 



W<p = a<p — — m u N. 

n M p i 

We have included a numerical factor oca, in the definition of Wa, given above. Although in the uniform density 
approximation we have ola, = 1, we chosen to leave ola, in the above equation so that our results can be more easily 
reassessed in light of the more accurate evaluation of Wa, performed in appendix [Dj The total energy is then: 



E(R) 



N ((N 2 / 3 /R 2 + mfj V2 - Nra^j 
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and p oc 1/R 3 . We will have equilibrium when dE(R) / dR — which implies: 



JVl/3 



(TV 2 / 3 + m 2 e R 2 ) 
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3m 2 u N 2 / 3 
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\n\ 



t RN~ 1/3 . 



We note that the term on the left-hand side of the above expression is always less than 1, and that both terms on the 
right hand side are positive definite. This implies that there is a maximal value of N. This maximal value is found by 
setting the left hand side to 1 and solving for N. Because both terms on the right hand side are positive definite, the 
maximum value of N, with a chameleon field present, will be less than or equal to the value it takes in pure General 
Relativity . We can see that, both with an without a chameleon, we have: 



N < N max 



5^, 
3m?, 



3/2 



Following what was done for the braneworld corrections to gravity in [58| we define: 

R = 7V 1/3 /™ e , 

and x = R/Rq. We also define po = 3m u N / '4ttRq = 3m u m\ / An , and Y — (N/N ma 
dE/dR = 0, is therefore equivalent to: 
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Hydrostatic equilibrium, 
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This is the mass-radius relationship for a white-dwarf star. The Chandrasekar mass- limit follows from setting x — 0. 
We can see that it is unchanged by the presence of a chameleon field. The chameleon field will, however, alter the 
mass-radius relationship. We shall assume, and later require, that the second term on the right hand side, i.e. the 
chameleon correction, is small compared to the first. Solving perturbatively under this assumption we find: 
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The maximum value of x for relativistic white-dwarfs (j>f > m e ) is x = 1. In order for our assumption that the effect 
of the chameleon field could be treated perturbatively to be valid, we need: 
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The right hand side is clearly increasing with x, and so we evaluate it at x = 1. For corrections to the mass-radius 
relationship to be smaller than 10% we must therefore require: 



i\/2aA 
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Evaluating this expression, for a white dwarf, we find 

'3.19/3 x 10 8 



5.065 x 10 



< 0.1. 



M 
1 mur 
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i+i 



A^+t < o.l, 



The above expression provides us with an upper-bound on [3. For natural values of M and A, this upper-bound is 
strongest for n = —4 theories. 

For the corrections to the mass-radius relationship to be smaller than 10% when n = —4 we need 

f3 < 6.70A~ 1/4 x 10 18 . 

Alternatively, for the corrections to be smaller than 1% we require (3 < 1.19A -1 / 4 x 10 18 . As the data improves it 
might, in future, to be able to limit any such corrections to being smaller than 0.1%. In this case we would need 
(3 < 2.18A~ 1/4 x 10 17 . In evaluating these limits, we have used the accurate value for a c f,(n — —4) found in appendix 
ID] acf,(n = —4) = 0.58. Despite the fact that these represent some of the best upper bounds on /3 for chameleon 
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fields, it is clear that j3 1 is still allowed. Even if the astronomical data improves to the point where we can rule 
out corrections at the 0.1% level, we would still be unable to rule out M p ;//3 > ITeV. 

The calculation for a neutron star proceeds along similar lines. In a neutron star, the neutrons provide both the 
degeneracy pressure and the gravitational potential. We must therefore replace both m u and m e by m n . For the 
correction for the mass-radius relationship to be less than 10% we must require: 

(8.15 x 10" 51 ) a* ^ — = r A— <0.1. 

v \ \n\ J \1 mm -1 / 

71+4 

The left hand side of the above expression is a factor of (m e /w n )™+ 1 smaller than the equivalent expression for a 
white dwarf. It follows that, for all n = —4, a weaker bound on the parameters results. When n = —4 we find the 
same bound. In FIG. [5] we have plotted the white-dwarf bounds on 0, M and A. We have, conservatively, assumed 
that corrections to the mass radius relationship are smaller than 5%. The plots are similar to those done for the 
Eot-Wash and WEP bounds. The whole of the shaded region is allowed and the plots for other theories with n < —4 
are similar to those with n = —8. Similarly, the plots for other n > theories looks much the same as the n — 4 plot 
does. 

These white dwarf bounds are included in the plots of section HXl where all the bounds on these chameleon theories 
are collated. 



B. General Relativistic Stability 



We have already derived conditions for the effect of the chameleon to be small compared to that of the Newtonian 
potential and thus produce only a negligible change to the mass-radius relation. It also is necessary to consider how 
the inclusion of a chameleon affects the stability of white dwarfs and neutron stars. This requires the inclusion of 
general relativistic effects. 

A compact body in hydrostatic equilibrium (d_E/di? = 0) will be stable against small perturbations whenever 
d 2 E/dR 2 > 0, i.e. we are at a minima of the energy. For a proof of this result see [13, H^. The onset of instability 
occurs when d 2 E/dR 2 = 0. For Newtonian gravity this occurs when the star becomes relativistic i.e. T — 4/3. 
Ignoring general relativistic effects but including chameleon effects we find that d 2 E/dR 2 is given by: 



d 2 E _ 2dE N 4 ' 3 f N l / 3 m 2 e R 2 /?i + 4\3a 
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When dE/dR — 0, the contribution from the chameleon field to the right hand side of this equation is positive. It 
follows that the effect of the chameleon field is to increase the stability of white dwarfs and neutron stars i.e. it makes 
d 2 E/dR 2 more positive. 

It is well known, jsjj], that General Relativity alters the stability of white dwarfs and neutrons stars. When GR 
effects arc included gravity is generally stronger. As a result, it tends to have a destabilising effect on configurations 
that are stable when studied in Newtonian physics. In the absence of chameleon corrections, but including GR effects 
(assuming GM s t ar /R <§; 1) the criterion for stability is roughly: 

where k ~ 0(1), [57[. 

Whilst GR destabilises white-dwarf stars, we have just noted that the chameleon force acts to stabilise them. In 
this section we will see how the leading order general relativistic effects balance out against the chameleon force, and 
study their cumulative effect on the stability of compact objects. The assumption that general relativistic effects 
are small means that these results will be more accurate for white dwarfs than they will for neutron stars. A full 
derivation of the potential energies associated with the leading order general relativistic effects can be found in [57l ]. 
For the sake of brevity we shall not repeat that analysis here but merely quote the results. 

We shall assume that the electrons in the white dwarfs are approximately relativistic and so satisfy P — Kp 4 ' 3 . 
We shall also assume spherical symmetry. Just as we do in appendix [Dl we evaluate the different contributions to 
the energy of the white-dwarf under the assumption that the fluid satisfies the Newtonian equation of hydrostatic 
equilibrium at leading order. We then solve the resulting Lane-Emden equation for P numerically. This procedure is 
valid if one only wishes, as we do, to calculate the GR and chameleon field corrections to leading order. 



41 



The internal energy of a n = 3 polytrope is: 

U = hKpl^muN, 

and the gravitational potential of the star is given by: 

W = -k 2 p 1 J 3 M- 2 (m u N) 5 ^, 

where k\ and k 2 are found in ref. [57} to be (for T = 4/3): 

fci = 1.75579, k 2 = 0.639001. 

In addition to the correction coming from General Relativity (which we will consider shortly) we also need to account 
for the fact that the electrons are not entirely relativistic. At leading order, this gives the following correction to the 
internal energy: 

AU = hml/ifiem^muNp- 1 / 3 , 
. Finally, the leading order general relativistic contribution to the energy is found to be of 



where fc 3 = 0.519723, 
the form: 



AW G r = 



-k 4 M-\m u N)y 3 n^ 



where U = 0.918294, 151 



'GR = -K4ivi pl ym u i\ ) ■ p c 
Including the effect of the chameleon, the energy of the white dwarf is given by: 
E = U + W + AU + AW G r. + W<f,. 

We define B = faK, C = k 2 M pl 2 , D = k 3 m 2 Jm 2 J 3 and F = k 4 M~^. We also define H = a <t> /\n\\f3(t){p c )M- l 1 \p 1 J {n+1) . 
As defined, H is actually independent of p c . With these definitions, the energy is extremised when: 
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At leading order, we drop the terms proportional to D, F and H and recover the standard Chandrasekar limit for 
the mass of a white dwarf: 



M stB 
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M, 
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where Mq is the mass of the Sun. Instability begins to occur when d 2 E/d 2 p c 
equivalent to: 



which, given that dE/dp c — 0, is 



o V Pc ~o' tM starPc H ~~T H Pc 

3 3 n + 1 



1+1 = 0. 



(42) 



Solving this equation for p c gives a critical density, p cr n, at which instability occurs. We find that, for all n =/= —4, 
there will be a chameleon induced correction to p cr it- When n ^ —4, we must either have n < —6 or n > and so then 
(n + 2)/(n + 1) > 4/5 > 1/3. This observation means that eq. |42|) still has solutions. The effect of the chameleon is 
to raise the value of p cr it- Even in pure General Relativity it turns out that this critical density is so high that it will 
only be important for white dwarfs in which the core is 4 He, [l?) • In all other cases, except for that where the core is 
12 C, Pcrit is greater than the ncutronisation threshold, and so such high-density white dwarfs will not occur. In the 
case of carbon white dwarfs, the central density is in fact limited by pyconuclear reactions, [571 I. Since the addition of 
a chameleon field raises p C riu this change can only be potentially important for Helium white dwarfs. 

If the effect of the chameleon is small then, at leading order, p cr n = CB 2 /DA 2 = 2.65 x 10 10 gem -3 for 4 He stars. 
The leading order chameleon correction to p cr it is 



Sp c 



3(n + 4) / a (2m u ) 2 / 3 ^ 



Pcrit 




(43) 



If we wish to require that the presence of a chameleon do little to alter the stability properties of white dwarfs in 
general relativity, we will need Sp cr n/ Pcrit "C 1. This gives us another upper bound on j3. In general, however, it is 
not competitive with the white-dwarf mass-radius relation bound on f3. The requirement Sp cr it/ Pcrit < 0.1 is included 
in the plots of FIG. E 
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FIG. 8: [Colour Online] Constraints on chameleon theories coming from white dwarfs and neutron stars. The shaded area shows 
the regions of parameter space that are allowed assuming that any alterations to the white-dwarf mass-radius relationship are 
at the 5% level or smaller and that the chameleon only induces changes that are smaller than 10% in the maximum white-dwarf 
density, p cr it- Neutron star bounds are not competitive with the white-dwarf constraints. The solid black lines indicate the 
cases where M and A take 'natural values'. The dotted-black line indicates when M = := M p i/f3 i.e. when the mass scale 
of the potential is the same as that of the matter coupling. Other n < —4 theories are similar to the n = —8 case, whilst the 
n = 4 plot is typical of what is allowed for n > theories. The amount of allowed parameter space increases with \n\. 



C. Discussion 



It should be noted that, in this section, the bounds that have been derived on /3 have been found under the 
assumption that the chameleon field couples to relativistic matter in the same way as it does to normal matter i.e. it 
just couples to the rest mass energy density of matter. As we noted in the introduction, however, it is usually the case 
that the chameleon in fact couples to some linear combination of the energy density and the pressure e.g. p + ujP. In 
the simplest models ui = —3. The total energy density, ptot, of the star with equation of state P = KpQ is given by: 

P 

Ptot = Po + jT— j" = Pa+pP 

where pa is the rest mass energy density and T = 1 + 1/p. In calculations presented above, we have implicitly assumed 
that ptot + wP = po + (p + ui)P ~ pq. The bounds that we have derived come from the sector where the matter in the 
star is relativistic i.e. p = 3. If the chameleon couples to matter through the trace of the energy momentum tensor 
i.e. oj = — 3 then we do, in fact, have ptot = Po, just as we have assumed. 

In the non-relativistic case, P <C po, and so ptot + ojP ~ po is always true. Even in the relativistic case, since 
P/p = Kp\j Z <C 1 for white dwarfs, and P/p — Kp^ 3 ~ O(l) for neutron stars, we always have P < O(p), and so 
different values of u) will only alter our bounds by at most an 0(1) quantity. There is one caveat: if oj is chosen so that 
ptot + can become negative, then the n > chameleon field theories will cease to display chameleon behaviour. 
This would immediately rule them out for all /3 > 1. 
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VIII. COSMOLOGICAL AND OTHER ASTROPHYSICAL BOUNDS 
A. Nucleosynthesis and the Cosmic Microwave Background 



The compact object bounds present above constrain a chameleon field behaves in very high density backgrounds 
whereas cosmological bounds on chameleon theories constrain how the behaviour of the chameleon field in low-density 
backgrounds. We have assumed that the chameleon couples to the energy density and pressure of matter in the 
combination: 



In the radiation era P w 3p. Provided then that p(l + w/3) > 0, i.e. ui > —3, and (3 is large enough, the chameleon 
will simply stay at the minima of its effective potential, which is itself slowly evolving over time. For this to be the 
case it is required that: 
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where d> r is the value of d> at its effective minima: 
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where we have used H 2 — p / ZM^ and p — —AHp as appropriate for the radiation era. 

We shall show below that we must require that \(3(j) c /M p i\ < 0.1 since the epoch of nucleosynthesis. When (3 > 1, it 
follows from eq. (|4"5| that this requirement is enough to ensure that \Ocj) c \ <C /3p(l + uj/3)/M p i provided that u> > —3. 

The simplest, and perhaps the most natural way, for the chameleon to interact with matter in a relativistically 
invariant fashion, however, is for it to couple to the trace of the energy momentum tensor i.e. oj = —3. When uj = —3 
the above analysis does not apply. The strongest bounds on the parameters of chameleon theories arise in u) = —3 
case. 

When oj = —3 we must evaluate p — 3P. Although p « 3P in the radiation era, p — 3P is not identically zero. 
Following [18|, we find, for each particle species i: 

45 H 2 M 2 pl 9l 

3Pi = — „ 77^\ T \ m i/T), 



TT A 8tt g,{T) 



where g, = £ boS o„ s 5l boso "(T 4 /T) 4 + (7/8) E/ erm40 „ s g( ermlon {T t / Tf is the standard expression for the total number 
of relativistic degrees of freedom; g t and Tj are respectively the degrees of freedom and temperature of the iih particle 
species. The function r is defined by: 



t(x) = x 



du 



e u ± 1 



with the + sign for fermions and the — for bosons. This function goes like x 1 when i<1 and e~ x when x 3> 1, but 
it is 0(1) when x ~ 0(1). The case (3 ~ 0(1) and n > was analysed very thoroughly in 18]]. We will therefore 
restrict ourselves to looking at the n < — 4 cases. We will also consider, for all n, what new features emerge when we 
take f3 ^ 1. In both cases we will see that theories with n < —4 and/or (3 3> 1 are much less susceptible to different 
initial conditions than those with n > and (3 ~ 0(1). 

We consider the cases n < —4 and n > separately below. Before we do so, we note some similarities between the 
two cases. Whatever the sign of n, the effective potential will have a minima at: 



,(T) = M 



M p iM 3 n\ 
(3p(T) 
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where we define: 

p{t) = j> 4 - m . 

i 

From the form of pi — 3Pj, it is clear that p will be dominated by the heaviest particle species that satisfies mi < T. 
When rrii < T, r a (m^/T) 2 and so the largest value of rm dominates, whereas the contribution from species with 
rrii 3> T are exponentially suppressed. The function T(m,/T) is peaked near rrii/T = 1. 

If the chameleon is not at <f> m i n , and (4> m in/4>) n+1 <C 1, then this peak will result in the value of (f> being 'kicked' 
towards (f> m in{T). We label the distance that cfr moves due to this 'kick' by {A<p)i where 

(A*), „ - 



*( m ») 

The 7/8 is for fermions and the I for bosons. This formula is valid so long as \<t> m in(T) — <j)\ > \A(f>i\ i.e. so long as 
the kick is not large enough to move <j> to its minimum. 

The largest jump of this sort will occur for the smallest value of g%/g*. It will therefore occur when electrons 
decouple from equilibrium at T w 0.5 MeV. If, however, \<f> — 4> m in(T)\ is smaller in magnitude than this above 
quantity, then (A<p)i — 4> m i n — (f>. Whether or not (f> will stay at <t> m in{T), as it evolves with time, will depend on the 
mass of the chameleon at <fr m in(T). If m 2 ^ 3> H 2 then it will stay at the minimum. Otherwise it will tend to slowly 
evolve towards values of <p for which (4> m in / 4>) n+1 < 1- 

If {<Pmin/ 4>) n+1 > 1 t nen <p will very quickly (in under one Hubble time) roll down the potential. It will either 
overshoot (j) m in, or if the mass of the chameleon at (f> m in is large enough, stick at <j> m in- We can therefore assume that 
our initial conditions are such that, in the far past, <fi is either at 4> m in or (<fi m in / 4>) n+1 < 1- 

1. Case n < —4 

We note that when n < —4, before any jump, we have (\> m inl4> > 1 and so \(j) m in — <j>\ < \<j>min\- It follows that: 

^mi«(r«mi) _ 45(n+l)ff 2 (T^m t ) 
{A<p)i 7r 4 m 2 

(0mm) 

This seems to suggest that, if (A<f>)i is large enough to move <j> to 4> m in{T), then we will necessarily have m^epmin) > 
H 2 , and so <j> will stay at 4> m in(T) in the subsequent evolution. However, this is not quite the case. As T drops 
below rrii, the ith species decouples and its energy density decreases exponentially. This causes p, and consequently 
m<p(4>min(T)) , to decrease quickly. Roughly, p shortly after decoupling will be a factor of (raj/mi) 2 smaller than it 
was before, where the jth species is the most massive species of particle obeying rrij < rrii. 

If (f> reaches mm (T w mi) with the ith kick, then as T decreases the chameleon will roll quickly down to the 
potential towards mm (T < mi); |0 mm (T < mi) / <j) m in(T w m»)| < 1. For (f> to stick at 4> m in, for m <T < rrn, and 
not overshoot it, we need: 



H 2 (T) (T\~^ T (m i X^ 1 H 2 ( mi ) 



m|(0 mm (T)) \ rrii ) \m 3 ) rnr^ m in{mi)) 
Since T < mi in this region, and nj (n + 1) > 0, it is sufficient to require: 
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We know that m^ oc /3("+ 2 )/(™+ 1 ) and so the above condition is more likely to be satisfied for larger values of (3 than 
for smaller ones. 

We note that, even when <p ^ 4> m im we cannot have 4>/4> m in 3> 1. If this were the case initially, when <f> = <pi say, 
then the gradient in the (p potential would be very steep and in one Hubble time <fi would move a distance A<p where: 

A(j) a ( m i\ 2 ( M PI \ ( 0mi" X ( ' 
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It is clear from this expression that for (3 > 0(1), 4>/<fi m in will decrease very quickly and pretty soon 4>/<t>min < 1. 

When (3 > 0(1) it is therefore valid to assume that, for almost of all of the radiation era evolution, cj)(T) / cj) min (T) < 
1. The larger (3 becomes, the greater the extent to which this assumption holds true. 

The main purpose of the above discussion is to illustrate that, for theories with n < —4, the late time behaviour 
of <p will be virtually, independent of one's choice of initial conditions. The same is not true, or at least not true 
to the same extent, of theories in which n > 0. The corollary of this result is that \j3(j>/M p i\ could be very large 
at the beginning of the radiation epoch and still be less than 0.1 from the epoch of nucleosynthesis onwards. This 
would correspond to the masses of particles during the early radiation era being very different from the values they 
have taken since the epoch of nucleosynthesis until the present day. The larger (3 is, the easier it is to support large 
changes in the particle masses. This is one reason why (3 ^> 1 is theoretically very interesting. In chameleon theories 
with (3 1, the constants of nature that describe the physics of the very early universe (i.e. pre- nucleosynthesis) 
could take very different values from the ones that they do today. This could have some interesting implications for 
baryogenesis at the electroweak scale. 

The best radiation-era bounds on (3, M and A come from big-bang nucleosynthesis (BBN) and the isotropy of the 
cosmic microwave background (CMB). As noted in [HI], a variation in the value of <fi between now and the epochs of 
BBN and recombination would result in a variation of the particle masses (relative to the Planck mass) of about: 



Am 



M p i 

BBN constrains the particle masses at nucleosynthesis to be within about 10% of their current values, [Hj]. This 
requires: 

\<Pbbn\ < O.iP^Mpt. 

We have argued above that 4>bbn I <j>min{TBBN) ^ 1 and so the above condition will be satisfied whenever: 

\&mn{T B BN)\ < O.l/T^Mpl. 

Nucleosynthesis occurs at temperatures between 0.1 MeV and 1.3 MeV. Since 4> m in(T) decreases with temperature, 
we conservatively evaluate the above condition with Tbbn = 2 MeV. At such temperatures the largest contribution 
to p will come from the electrons (with m e (today) = 0.511 MeV) and: 

. _ g e T 2 ml _ ge (MeV) 4 
P ~ 24 24 

where g e = 4 (2 from the electrons and 2 from the positrons). The BBN constraint on f3, A and M for n < —4 is 
therefore: 

. { 1.8/3 x 10 6 \^ / m 

An — rrn — rr <Ux 10 37 . 

\ A|n| / \ 1 mm 1 / 

This is, however, a weaker bound on {(3, M, A} than the white dwarf mass-radius relation constraint discussed in 
section I VIII above. 

Another important restriction on these chameleon theories comes out of considering the isotropy of the CMB 49] . 
As is mentioned in 18] a difference between the value of (f> today and the value it had during the epoch of recombination 
would mean that the electron mass at that epoch differed from its present value Am e /m e m (3A(j)/M p i. Such a change 
in rrc e would, in turn, alter the redshift at which recombination occurred, z rec : 
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WMAP bounds z rec to be within 10% of the value that has been calculated using the present day value of m e , \18\ - 
We define p rec and pbbn to be, respectively, the value that p takes at the recombination and BBN epochs. Now 
Prec ^ Ptoday, where ptoday is the cosmological energy density, and </> cx / 9 _1 /( n + 1 ) j it follows that \A<f>\ s» \4> r ec\ for 
rc < —4. <}> rec is the value of <p had during the epoch of recombination. However, since p rec < pbbn, we also have 
<j>rec/ 4>bbn < 1, and so this CMB bound is always weaker than the one coming from BBN. For this reason we do not 
evaluate the CMB constraint here. 
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2. Case n > 



When we considered theories with n < —4, the analysis was made easier by the fact that V(<f>) had a minimum. 
This allowed us to bound the magnitude of <fi to be approximately less than that of 4> m in(T). However, when n > 0, 
the potential is of runaway form and we cannot bound the ma gnit ude of (f> in such a way. Many of the issues involved 
with n > potentials were discussed very thoroughly in ref. In that work it was assumed that (3 ~ 1 and so, 

generically m^,(4> m in) *C H in the radiation era, due the fact that p -C p. However, if (3 <C 1, it is possible for to 
be large compared to H. 

As with the n < —4 case, it is not necessary to require that <f> is at its effective minima during the whole of 
the radiation era. All that is really needed is for ft be sufficiently close to the minima at the epochs of BBN and 
recombination that we are able to satisfy their constraints. It was shown in 18] that the total sum of all of the kicks 
that occur before BBN will move the chameleon a distance of approximately: (3M p i. BBN requires that: 

|A0 BSA r| = fttoday ~ 4>BBN < 0.1/3 _1 M p ;. 

Provided that, at the beginning of the radiation era, ft <C 0M p i/8ir, then ft at BBN will easily satisfy the above bound 
provided that (3(j>toda V /M p i < 0.1. 

The first of these requirements is a statement about initial conditions. It is clear that the larger (3 becomes, the 
less restrictive this condition is. Indeed for (3 ^> 1, it is quite possible to have had (3ft/M p i ~ 0.1O(/3 2 ) S> 1 at the 
beginning of the radiation era and still satisfy this bound. The larger the matter coupling is, the less important the 
initial conditions become, and the greater the scope for large changes in the values of the particle masses and other 
constants to have occurred between the pre-BBN universe and today. 

Recombination enforces a similar bound to BBN: 

A^-rec ^ Pifttoday ftrec) < Q ^ 
Zrec Mpl 

If the requirements on the initial conditions are satisfied, we will have fttoday 3> ftrec, 4>bbn and so, in both cases, 
Aft w fttoday- We must therefore require that 

Pfttoday/M pl < 0.1. 

fttoday oc p c 1 ^ n+1 * ) where p c is that part of the cosmological energy density of matter that the chameleon couples 
to. For the most part, we have, in this paper, assumed that the chameleon couples to all forms of matter with the 
same strength. However, up to now, we have only been concerned with baryonic matter. Even if the chameleon to 
baryon coupling is virtually universal, there is not necessarily any reason to expect the chameleon to couple to dark 
matter with the same strength. It is possible that the chameleon does not interact with dark matter at all. Since it 
clear that fttoday is a very important quantity when it comes to bounding n > chameleon theories, it is crucial to 
know what fraction of the cosmological energy density the chameleon actually couples to. The smaller p c is, the larger 
4>today will be. The larger the value of fttoday , the more restrictive the resultant bound on {(3,M, A}. This implies 
that cosmological bounds on a chameleon theory that couples only to baryonic matter will be stronger than those on 
a chameleon that also couples to dark matter. Not knowing how the chameleon couples to dark matter, we chose to 
be cautious, and err on the side of specifying a bound that is perhaps slightly too tight, rather than too loose. We 
therefore work on the assumption that the chameleon only couples to baryonic matter, and so p c ~ 0.42 x 10 - 30 g cm" 3 , 
|49j . Under this assumption, we find: 
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8.34 x 10- 3 M I - J I — — T I < 1. 
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If the chameleon to dark matter coupling is similar in magnitude to the baryon coupling then p c i=s 2.54 x 
10 -30 g cm -3 and we have the less restrictive bound: 



1.17 x IQ-^I3\ " +1 / M 
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The bound that we have just found has come about from the requirement that the particle masses at BBN and 
recombination are within 10% of the values they take in regions of cosmological density i.e. ptoday ~ 10~ 30 gcm~ 3 . 
However, all cosmological determinations of the particle masses, and indeed of the other constants of nature, have 
come from analysing measurements made in regions with densities much greater than the cosmological one. For 
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instance, recent cosmological determinations of m p /m e have employed the absorption and emission spectra of dust 
clouds around QSOs [6Ul67|. These dust clouds have typical densities of the order of p ~ 1CT 25 — 10~ 24 gcm~ 3 . If 
we take p c to be 10 -25 gem -3 , then the BBN and CMB bounds would only require: 

4.59 x l(T 24 /?\ ^ / M \^ 
n J \lrnm _1 y 



3.52 x 10" 



3. Summary 



We have plotted the BBN and CMB constraints on {/3, MX} in FIG. M We have plotted what occurs in the most 
restrictive case i.e. when the chameleon couples only to baryons. The whole of the shaded region is currently allowed. 





FIG. 9: [Colour Online] Constraints on chameleon theories coming from the particle masses at big-bang nucleosynthesis and 
the constraints the redshift of recombination. The shaded area shows the regions of parameter space that are allowed by the 
current data. The solid black lines indicate the cases where M and A take 'natural values'. The dotted-black line indicates 
when M = M<f> := M p i//3 i.e. when the mass scale of the potential is the same as that of the matter coupling. Other n < —4 
theories are similar to the n — —8 case, whilst the n — 4 plot is typical of what is allowed for n > theories. The amount 
of allowed parameter space increases with |n|. In these plots we have assumed that the chameleon couples only to baryons. 
Slightly weaker constraints result if the chameleon additionally couples to dark matter. 

Another class of potentially important cosmological bounds can be derived by employing astronomical bounds on 
the allowed variation of the fundamental constants of nature during the matter era. We discuss this further below. 



B. Variation of fundamental constants 



Any variation in a chameleon held will lead to a variation in the masses of the particle species to which the chameleon 
couples. This variation is relative to the fixed Planck mass, M p i = G~ 1//2 . If the chameleon couples to all matter 
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particles in the same way, then all the fundamental particle masses will vary in the same fashion and so their ratios 
will remain constant. It is also feasible to construct theories where a variation in the chameleon leads to a variation 
in some other fundamental 'constants' of nature. For instance one might propose a theory where the fine-structure 
constant, a em , is given by a em = f(f3(f>/M p i) for some function /. If this is the case then, in addition to bounds on 
any allowed variation in the particle masses, we would also have to apply the whole range of very stringent bounds 
on the possible variation of a em mentioned above. 

For the purposes of this section, however, we assume that a em is a true constant. It should be noted that even 
if a em does vary with <f> at the same level as the particle masses do, the resulting bounds on the parameters of the 
theory are only slightly more stringent than those already found. 

The best matter era bounds on the variation of the particle masses come from measurements of the ratio /i = m p /m e , 
(l3l l6ll . We shall assume that the chameleon couples to protons with strength (3 P and to electrons with strength 
(3 e . The relative change in the proton and electron masses is then given by: 

Am p f3 p A(j) 

m p M p i 
Am e _ /3 e Ac/) 

m e M p i 

Without any a priori knowledge about the magnitude, or sign, of f3 e — (3 P it is difficult to derive any bounds on 
chameleon theories simply by considering A/1//1, where fi = m p /m e . The simplest assumption one could make about 
the matter coupling, f3, is that it is universal i.e. j3 p = (3 e = /?■ If this is the case then Afi — identically. An 
alternative, but still very reasonable, assumption about the chameleon coupling, which would produce Afj, ^ 0, is 
that the chameleon couples to all fundamental particles with the same universal coupling, flu, but that the QCD 
scale, Aqcd, is independent of 4>. When the quark masses vanish the proton mass is proportional to Aqcd- The 
masses of the three lightest quarks, m Ul and m s , are considerably smaller than Aqcd and as a result contribute 
only a small correction to the proton mass (at approximately the 10% level). If Aqcd is (^-independent, we expect 
the proton mass to depend only weakly on 0, through the quark masses, and so fl p ~ 0(flu /10) whereas fl e = flu- 
Since the mass of the electron is so much smaller than the proton mass, the coupling of the chameleon to baryons is 
approximately given by fl p , and so it is fl p that is constrained by the experiments mentioned in sections I VII and IVIII 
However, the BBN and CMB requirements constrain fl e — flu- A change in </> of Acj) would therefore induce a change 
in p of: 

*a « -9^. (45) 

p M p i 

As we reported above, the Reinhold et ai, [6l| . result is consistent with a difference between the laboratory value of 
p and the one measured in such a dust cloud at the 3.5a level: Ap/p = 2.0 ± 0.6 x 10 -5 where Ap = pdust — piab- 
It should be noted that in the context of the chameleon theories considered here A<f> = <f>dust — 4>iab is always positive 
and so A/i//i < if j3 p = fl e /10 as we have assumed. Under this assumption, it is not possible to reproduce the result 
of Reinhold et al. If we had alternatively assumed that the fundamental particle masses are true constants but that 
Aqcd is ^-dependent, then we would be able to have Ap/p > 0. 

We interpret the Reinhold result conservatively i.e. as limiting any variation in p to be at or below the 2 x 10~ D 
level. We shall also assume that the laboratory experiments that measure p are performed close enough to the Earth's 
surface that the background value of <p for these experiments is approximately the value the chameleon takes inside 
the Earth. This is last assumption is also a conservative one i.e. it will result in a tighter bound on the chameleon 
theory parameters. Taking the density of Earth to be p Earth ~ 5.5 gem -3 and the density of the dust clouds from 
which the absorption spectra come to be pdust ~ 10~ 25 g cm -3 we find that we must require: 

2.88 x 10-M ( 2 ^f) * (r^j) m X^<l, n<-A 

,, , /4.59 x 10- 2i B\ *r / m \ ^ 

1.58 x 10" 3 n -) —) <1, n>0. 
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For fixed M and A this places an upper-bound on 0. When n > 0, the bounds coming from varying-/^ are competitive 
with the other cosmological constraints and they provide a weak upper bound on (3. When n < —4, however, the 
white dwarf bounds of section IVIII still provide the best upper bound on (3. It should be noted that the bounds on 
{/?, M, A} deduced from measurements of A/i//i are highly model-dependent. For this reason we do not plot them 
here. 



\ |n| / \ 1 mm 
(3 < 1.42A~ 1/4 x 10 19 , 
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IX. COMBINED BOUNDS ON CHAMELEON THEORIES 

All of the chameleon theories considered in this work have a two-dimensional parameter space, spanned either by 
M and (3 (n ^ —4), or by A and (3 (n = —4). We combine the constraints found in sections IVHIVIIII to bound the 
values of (3 and M (or A) for different n. We plotted the constraints for n — —8, —6, —4, 1 4 and n — 6 in FIG. [TO] In 
these figures we have included all the bounds coming from the Eot-Wash experiment, [H, [3£j, as well as those coming 
from Casimir force searches, [4(], We also include the bounds (labeled Irvine) coming from another search for 
Yukawa forces, [65]. We also show how current WEP violation experiments constrain these theories - i.e. experiments 
that have actually been done as opposed to the putative WEP violation experiment described in section IVI CI The 
white-dwarf and BBN constraints are also included in the plots, however, for these theories, they are always weaker 
than those laboratory tests (when (3 < 10 20 ). The plots for other n < —4 theories are very similar to the n = — 8 
and n = —6 plots, whilst the n = 4 and n — 6 cases are representative of n > theories. In general, the larger |n| 
is, the larger the region of allowed parameter space. This is because, in a fixed density background, the chameleon 
mass, m c , scales as |n + l| 1 / 2 /|n| 1 / 2 ( n+1 ), and so m c increases with \n\. The larger m c is, in a given background, the 
stronger the chameleon mechanism, and a stronger chameleon mechanism tends to lead to looser constraints. 

We have indicated on each plot how the variety of different bounds, considered above, combine to constrain the 
theory. In each plot, the whole of the shaded area indicates the allowed values of M and f3 (or A and (3). Three 
satellite experiments (SEE [53|, STEP and GG (55[) are currently in the proposal stage, whilst a fourth one 
(MICROSCOPE [Hi]) will be launched in 2007. These experiments will be able to detect WEP violations down to 
i] = 10~ 18 . The more lightly shaded region on the plots indicates those regions of parameter space that could be 
detected by these experiments - we have assumed that their WEP violation experiments to be along the lines of the 
putative experiment described in section fVI CI 

In n > theories, the increased precision, promised by these satellite tests, is not the only advantage that they offer 
over their lab-based counterparts. In space, the background density is much lower than in the laboratory. As a result, 
the thin-shell condition, eq,. (|19b ). is generally more restrictive for bodies in the vacuum of space than it is for the 
same bodies here on Earth. It is therefore possible for test-bodies, that had thin-shells in the laboratory, to lose them 
when they are taken into space [13, EH, Hl| . If such a thing occurs for the test- masses used in the aforementioned 
satellite experiments, then it is possible that they might see WEP violations in space at a level that had previously 
been thought ruled out by Earth-based tests i.e. 77 > 10~ 13 . This interesting feature of n > chameleon theories was 
first noted in [25j . It is important to stress that this effect will not occur if (3 is so large that the satellites themselves 
develop thin-shells [13, [HI]. ^ n chameleon theories where the potential has a minimum, i.e. n < —4, the thin-shell 
condition, eqs. (TiT)k & c), is only weakly dependent on the background density of matter. As a result, n < —4 theories 
will generally be oblivious to the background in which the experimental tests of it are conducted. It is for this reason 
that future space-based tests will be better able to constrain n > theories than they will n < —4 ones. 

In all of the plots, we show (3 running from 10~ 10 to 10 20 . (3 < 10 -10 will remain invisible to even the best of 
the currently proposed space-based tests, and (3 > 10 20 corresponds to M p i/[3 < 500 MeV. The region in which 
MpiJ j3 < 200 GeV is, in fact, probably already ruled out. If (3 were so large that M p i/j3 < 200 GeV, then we would 
probably have already seen some trace of the chameleon in particle colliders. This said, without a quantum theory 
for the chameleon it is hard to say how chameleon theories behave at high energies. A result a detailed calculation of 
the chameleon's effect on scattering amplitudes is not possible at this stage. A full quantum mechanical treatment of 
the chameleon is very much beyond the scope of this work, but remains one possible area of future study. 

The chameleon mass (m c ), in a background of fixed density, scales as A"+ 1 M 2 (™ +1 > . As we mentioned above, 
the larger m c , is the easier it is to satisfy the thin-shell conditions, eqs. (fl9k -c). and the stronger the chameleon 
mechanism becomes. Since (n + 4)/(2(n + 1)) > and n/(n + 1) > 0, for all theories considered here, the chameleon 
mechanism becomes stronger as M — > 0, or A — > 00, and all of the constraints are more easily satisfied in these limits. 
It is for this reason that we truncate our plots for some small M and, when n = —4, for a large value of A. Values of 
M that are smaller than those shown, or values of A that are larger, are still allowed. 

The upper limit on M (and lower limit on A) has been chosen so as to show as much of the allowed parameter space 
as possible. When (3 is very small, the chameleon mechanism is so weak that, in all cases, <f> behaves like a standard 
(non-chameleon) scalar field. When this happens, the values of M and A become unimportant, and the bounds one 
finds are on (3 alone. This transition to non-chameleon behaviour can be seen to occur towards the far left of each of 
the plots. 

It is clear from FIG.HOlthat j3 3> 1 is, rather unexpectedly, very much allowed for a large class of chameleon theories. 
We can also see that, rather disappointingly, future space-based searches for WEP violation, or corrections to 1/r 2 
behaviour of Newton's law, will only have a small effect in limiting the magnitude of (3. If M p i/(3 ~ 1 TeV is feasible, 
pending a detailed calculation of chameleon scattering amplitudes, that chameleon particles might be produced at 
the LHC. 

The solid black line on each of the plots indicates the 'natural' values of M and A i.e. M ~ (0.1mm) -1 and 
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FIG. 10: [Colour Online] Combined constraints on chameleon theories. The whole of shaded area shows the regions of parameter 
space that are allowed by the current data. Future space-based tests could detect the more lightly shaded region. The solid 
black lines indicate the cases where M and A take 'natural values'. For n 7^ —4, a natural value for M is required if the 
chameleon is to be dark energy. The dotted- black line indicates when M = M<f, := M p i/f3 i.e. when the mass scale of the 
potential is the same as that of the matter coupling. Other n < —4 theories are similar to the n = —6 and n = —8 cases, whilst 
the n = 4 and n — 6 plots are typical of what is allowed for n > theories. The amount of allowed parameter space increases 
with \n\. 



A = 1/4!. The Eot-Wash experiment and measurements of the Casimir force rule out A = 1/4! in <fi 4 , however (expect 
when n = -6) M ~ (0.1mm)- 1 is allowed for all 10 4 < /3 < 10 18 . In particular, M pl //3 ~ 10 15 GeV w M GUT , i.e. the 
GUT scale, and M p i/(3 - 1 TeV are allowed. 

The dotted black line indicates the cases where M = M p i/j3 i.e. when there is only one mass scale associated with 
the chameleon theory. It is clear, however, that no such theories are allowed if /3 < 10 20 . In all cases we must require 
M <C M p i/(3. As noted in (l7l [25j], this requirement introduces a hierarchy problem in the chameleon theory itself. 
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This problem is not, however, present in 4 (n — —4) theory. 

The larger (3 is the stronger the chameleon mechanism becomes. A strong chameleon mechanism results in larger 
chameleon masses, and larger chameleon masses in turn result in weaker chameleon-mediated forces. A stronger 
chameleon mechanism also increases the likelihood of the test masses, used in the experiments consider in section fVTl 
having thin-shells. Large values of (3 cannot, therefore be detected, at present by the Eot-Wash and WEP tests, if A 
and M take natural values. It is for this reason that, as can be seen in FIG. [TQl the Eot-Wash experiment and the 
WEP violation searches place the greatest constraints on the parameter space in a region about (3 ~ 0{1) . Casimir 
force tests are much better able to detect large /3, but they can ultimately only place an upper-bound on M or 1/A. 

For large j3, Casimir force test succeed where the Eot-Wash experiment fails because in the former does not use a 
thin metal sheet to cancel electrostatic effects. If it were not for the presence of the BeCu sheet, the Eot-Wash test 
would be able to detect, or rule out, almost all [3^> 1 theories with M ~ (pa) 1 ^ 4 - 

Allowed couplings and force ranges for linear theories 
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FIG. 11: Allowed couplings and force ranges for scalar fields which couple to baryon number and have approximately linear 
field equations. The shaded area shows the parameter space that is allowed by the current bounds on scalar field theories with 
linear field equations. A = 1/m is the range of the force, and m is the mass of the scalar field. /3 is the matter coupling of the 
scalar. /3 > 1 is ruled out for all but the smallest ranges (currently A < 10~ 4 m). We refer the reader to [3^1 for a plot of the 
allowed parameter space for A < lCP 3 m. Irvine, Seattle and Boulder refers to |65|. [o2l. lo3l. loll] and [Hj] respectively. 

The fact that we can have (3 3> 1, in scalar theories with a chameleon mechanism, is entirely due to the non-linear 
nature of these theories. Almost all the quoted bounds on the coupling to matter are for scalar field theories with 
linear field equations. In such theories </> evolves according to: 

-U<j> = m c (j>+——, 
M p i 

where the field's mass, m, is constant (i.e. not density dependent). The <?!>-force between two bodies with masses Mi 
and M2, which are separated by a distance d, takes the Yukawa form 

_ _ (3 2 (l + d/X)e- d / x M 1 M 2 

■T12 — 



M p V 2 

where A = l/m c is the range of the force. The best limits of A and (3 come from WEP violation searches, [III H3, 
Ir32l . Ir33l . HH , and searches for corrections to the 1/r 2 behaviour of gravity, [38| . The 95% confidence limits on m and 
j3 for such a linear theory, where the field couples to baryon number, are plotted in FIG. [Tl] with the allowed regions 
shaded and the excluded regions left white. It is clear that (3 > 1 is ruled out for all but the smallest ranges (currently 
A < l(T 4 m = 0.1mm). 
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Cosmological force ranges for n = -8 theories Cosmological force ranges for n = -6 theories 
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Cosmological force ranges for n = 4 theories Cosmological force ranges for n = 6 theories 
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FIG. 12: [Colour Online] Allowed couplings and cosmological force ranges for chameleon theories. The shaded area shows the 
allowed parameter space with all current bounds. A cos is the range of the chameleon force in the cosmological background with 
density p ~ 10 -29 gcm -3 . It is related to the cosmological mass of the chameleon, m" s , by A cos = l/m^, os . The solid black 
lines indicate the cases where M and A take 'natural' values. Plots for theories with n < —4 or n > are similar to the cases 
n = —8, —6 and n = 4, 6 respectively. 

To make the comparison with the linear case more straightforward, we have replotted the allowed parameter space 
for chameleon theories with n = —8, —6, —4, 4 and n = 6 in terms of its coupling to matter, (3, and the range of 
the chameleon force cosmologically, A cos , and in the solar system, A so ;, in FIGS. [T2l and [T51 FIG. [T2l shows the 
cosmological range, whilst X so i is shown in FIG. [13] The solid black line, in each of these plots, indicates the case 
where M and A take their 'natural' values. We can clearly see that, in stark contrast to the linear case, chameleon 
theories are easily able to accommodate both /3 > 1 and A > lm. This underlines the extent to which non- linear 
scalar field theories are different from linear ones, and the important role that is played by the chameleon mechanism. 



53 



Solar system force ranges for n = -8 theories Solar system force ranges for n = -6 theories 
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Solar system force ranges for n = 4 theories Solar system force ranges for n = 6 theories 
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FIG. 13: [Colour Online] Allowed couplings and cosmological force ranges for chameleon theories. The shaded area shows the 
allowed parameter space with all current bounds. X so i is the range of the chameleon force in the solar system, where the average 
density of matter is p ~ 10 -24 gcm -3 . It is related to the cosmological mass of the chameleon, m s c ° l , by A so ; = l/m s c ° l . The 
solid black lines indicate the cases where M and A take 'natural' values. Plots for theories with n < —4 or n > are similar to 
the cases n = —8, —6 and n = 4, 6 respectively. 

X. CONCLUSIONS AND DISCUSSION 

In this article we have investigated scalar field theories in which the field is strongly coupled to matter. In particular, 
we have studied the so called chameleon scalar fields. A scenario presented by Khoury and Weltman [17] that employed 
self-interactions of the scalar-field to avoid the most restrictive of the current bounds on such fields and its coupling. 
In the models that they proposed, these fields would couple to matter with gravitational strength, in harmony with 
general expectations from string theory, whilst, at the same time, remaining very light on cosmological scales. In this 
work we went much further and show, contrary to most expectations, that the scenario presented in (l7| allows scalar 
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fields, which are very light on cosmological scales, to couple to matter much more strongly than gravity does, and yet 
still satisfy all of the current experimental and observational constraints. 

Previous investigations on such scenarios, [H, HH, IH, noted that an important feature of chameleon field 
theories is that they make unambiguous and testable predictions for near-future tests of gravity in space. This is 
timely as there are currently four satellite experiments either in the proposal stage or due to be launched shortly 
(SEE [H|, STEP [HI, GG [El] and MICROSCOPE [Hi]). A reasonably sized region of the parameter space of the 
chameleon theories considered here will be visible to these missions. Theories with very large couplings /5 > 1 will, 
however, remain undetectable. The ability of these planned missions to detect large (3 theories could, however, be 
exponentially increased if the experiments they carry were to be redesigned slightly in the light of our findings. 

Previous studies claimed that typical test masses in the above satellite experiments do not have a thin shell. 
Therefore, the extra force is comparable to their gravitational interaction. The chameleon model hence predicts that 
MICROSCOPE, STEP and GG could measure violations of the weak equivalence principle that are stronger than 
currently allowed by laboratory experiments. Furthermore, the SEE project could measure an effective Newton's 
constant different, by order of unity, from that measured on Earth. We have seen, in this paper, that both of these 
features are very much properties of chameleon theories with runaway (n > 0) potentials. They will not, in general, 
occur if the chameleon potential has a minimum (e.g. n < —4 theories). These features are also very much associated 
with a gravitational strength chameleon coupling i.e. (3 ~ 0(1), and will not occur if ft 3> 1, or j3 -C 1. 

The major result, presented in this work, is that current experiments do not limit the coupling of the chameleon to 
matter, (3, to be order 0(1) or smaller. Indeed, if we wish to have a 'natural' value of M in a V oc (p^ 1 theory then 
we must require (3 > 10 4 (or (3 < 10 -3 . If (3 3> 1, the test-masses in the planned satellite experiments will still have 
thin-shells. As such SEE, STEP, GG and MICROSCOPE, as they are currently proposed, will be unable to detect 
the chameleon and place an upper-bound on {3. 

We have shown that upper-bounds on the matter coupling, (3, can be derived from astrophysical and cosmological 
considerations. Also, if (3 is very large, of the order of 10 17 or greater, then it might even be possible to detect the 
effect on the chameleon on scattering amplitudes in particle colliders. This possibility is one avenue that requires 
further in depth study. 

We noted in the introduction that ^ > 1 could be seen as being pleasant in light of the hierarchy problem 
[l06l Il07l Il08j ] . If a chameleon with a large (3 where detected, it would imply that new physics emerges at a sub- 
Planckian energy-scale: = M p i/ (3. 

A large value of the matter coupling is also preferable to an order unity value in that it leads to the late time 
behaviour of the chameleon being much more weakly dependent on the initial conditions, than it would be if (3 < 0(1). 

The magnitude of (3(A<fi = <fii — 4>2)/M p i quantifies the relative amount by the particle masses differ between a 
region where <f> = 4>\ and one where <j) — 4>2- The larger the coupling is then, the easier it is for there to have been 
a very large difference between the current values of particle masses, and the values that they had in the very early 
universe (i.e. pre-BBN). If the particle masses were very different from their present values at, say, the epoch of the 
electroweak phase transition, then the predictions of electroweak baryogenesis could be significantly altered. 

In this paper have taken the chameleon potential to have a power-law form i.e V oc cf)~ n . This is certainly not the 
only class of chameleon potentials that it is possible to have. In general, any potential that satisfies (3V^ < 0, > 
and V^^/Vcf, > in a region near some (f> = <fro will produce a chameleon theory. In fact, a generic potential might 
contain many different regions in which chameleon behaviour is displayed. In some of these regions, the potential 
may appear to have a runaway form, and so behave qualitatively as an n > theory. In other regions, the potential 
might have a local minimum, leading to n < —4 type behaviour. The existence of the matter coupling provides one 
with a mechanism, along the lines of that considered by Damour and Polyakov [20| , by which the scalar field (f> can, 
during the radiation era, be moved into a region where it behaves like a chameleon field. The larger (3 is, the more 
effective this mechanism becomes. Given this mechanism, one important avenue for further study is to see precisely 
how general late-time chameleon behaviour is of a generic scalar field theory with a strong coupling to matter. 

In this paper we have avoided the temptation to linearise the chameleon field equation, eq. ([9]), when it is not 
valid to do so. We have, instead, combined matched asymptotic expansions with approximate analytical, and exact 
numerical, solution of the full non-linear equations to study the behaviour of chameleon field theories in more detail. 

The main results of this analysis were: 

• We found the conditions under which a body would have a thin-shell, and noted that the development of a 
thin-shell is related to the onset of non-linear behaviour. 

• We have shown that the far field of a body with a thin-shell is independent of the coupling strength (3: this is a 
generic property of all the chameleon field theories with a power-law potential. This /3-independence was seen, 
in section IVH to have important consequences for the design of experiments that search for WEP violations, and 
was seen to be vital in allowing theories with 3> 1 to be compatible with the current experimental bounds. 
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The (3 independence of the far-field of thin-shelled bodies was seen to cause the 0-force between two such bodies 
to be /3-independent also. 

• For (3 ~ 0(1) the best bounds on M and A currently come from the Eot-Wash experiment [H, [3§| For ^> 1, 
the best bounds on M and A were found to come from measurements of the Casimir force. 

• Non-linear effects were shown to limit the magnitude of the average chameleon mass in a thin-shelled body to 
be smaller than some critical value, m^ nt . Intriguingly, m^ rlt is independent of (3, M and A, and depend only 
on n and the microscopic properties of the thin-shelled body. 

• The experimental constraints on the coupling of chameleon fields to matter are much weaker than those on non- 
chameleon fields. In fact when n ^ —4, the constraints on large (3 theories are weaker than those in which the 
scalar field couples to matter with gravitational strength ((3 ~ 0(1)). Almost paradoxically, strongly-coupled 
scalar fields are actually harder to detect than weakly coupled ones. 

• Perhaps the most important result, though, is that the ability of table-top gravity tests to see strongly coupled, 
chameleon fields could be exponentially increased if certain features of their design could be adjusted in the 
appropriate manner. The detection, or exclusion, of chameleon fields with /3 > 1 represents a significant but 
ultimately, we believe, achievable challenge to experimentalists. Searches for large (3 chameleon fields represent 
one way in which table-top tests of gravity could be used to probe for new physics beyond the standard model. 
Whether, or not, chameleon fields actually exist, it is important to note those areas into which our current 
experiments cannot see and, if possible, design experiments to probe those areas. 

We have shown, in this paper, that scalar field theories that couple to matter much more strongly than gravity are 
not only viable but could well be detected by a number of future experiments provided they are properly designed 
to do so. This result opens up an altogether new window which might lead to a completely different view of the role 
played by light scalar fields in particle physics and cosmology. 
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APPENDIX A: PSEUDO-LINEAR REGIME FOR SINGLE-BODY PROBLEM 

In the pseudo-linear approximation, we assume that non-linear effects are, locally, everywhere sub-leading order. 
The cumulative, or integrated, effect of the non-linearities is not, necessarily small. This means that, whilst we assume 
that there always exists at least one self-consistent linearisation of the field equations about every point, we do not 
require there to be a linearisation that is everywhere valid. Instead we aim to contrast two linearisations of the field 
equations: the inner and the outer approximations to cj). 

The inner approximation is intended to be an asymptotic approximation to the chameleon that is valid both inside 
an isolated body, and close to the surface of that body. We take the isolated body to be spherically symmetric, with 
uniform density p c and radius R. Far from the body, r 3> R, the inner approximation will, in general, break down. 

The outer approximation is an asymptotic approximation to <f> that is valid for large values of r. We require that 
it remains valid as r — > oo. In general the outer approximation will not be valid for r ~ O(R). 

The boundary conditions on the evolution of <j) are: 

= 0. 

The first of these conditions is defined at r = 0. We will generally find that only the inner approximation is valid at 
r = 0. As a result we cannot apply the r = boundary condition to the outer-approximation. Similarly the r — oo 
boundary condition will be applicable to the outer-approximation but not to the inner one. 

Since we cannot directly apply all the boundary conditions to both approximations, there will generally be undefined 
constants of integration in both the inner and outer expansions for <f>. This ambiguity in both expansions can, however, 



dr 



= 0, 



d(j) 
dr 
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be lifted if there exists some intermediate range of values of r (r out < r < Ti n say) where both the inner and outer 
approximations are valid. 

Asymptotic expansions are locally unique (43l. El]. Thus, if both the outer and inner approximation are simultane- 
ously valid in some intermediate region, then they must equal to each other in that region. By appealing to this fact, 
we can match the inner and outer approximations in the intermediate region. In this way, we effectively apply all 
of the boundary conditions to both expansions. This method of matched asymptotic expansions is described in more 
detail in [45[ . 



1. Inner approximation 



Inside the body, < r < R, the chameleon obeys: 



d 2 (j, 2d<t> ^ r3 fM 3 \ n+1 /3 Pc 

' r^ = "^Af 3 — + ^-. (Al) 



dr 2 r dr \ <j> J M, 

The inner approximation is defined by the assumption 



M, 



Defining 



<\>c = M 

we see that the above assumption is equivalent to: 



i\M„,M 3 



n+l 



We define the inner approximation by solving eq. (|A1|) for <ft as an asymptotic expansion in the small parameter (5(0); 
we shall see below that S(r) < 5(0) := S. 

Whenever the inner approximation is valid we have: 

6M p i 

where the order 5 term is: 

n+l 



a ( r ) = *0£c r dr , r dr n T „ 



rM p i Jo J \<p(r") 

We have defined <f)(r) := <po + flpc^ 2 /QM p i for r < R. tfio is an undefined constant of integration. It will be found by 
the matching of the inner approximation to the outer one. For the inner approximation to remain valid inside the 
body we need: 

Outside the body, r > R, <j> obeys: 

f± + l^±_ _ nXM 3 (M^_\ n+1 
dr 2 r dr \ <p J 

Whenever S(r) < 1, we can solve the above equation in the inner approximation finding: 

~ 4>(r) + 4>s{r). 
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Outside the body, r > R, we define <j>(r) to be 

2M p i 3M p ir 

The order 5 term, 4>s(r), is given by the same expression as it was for r < R. 

The inner approximation will therefore be valid, both inside and outside the body, provided that: 

0(r) ^ 

In general, this requirement will only hold for r less than some finite value of r, r = ri n say. To fix the value of </>o, 
and properly evaluate the above condition, we must now consider the outer approximation. 



2. Outer Approximation 

When r is very large, we expect that the presence of the body should only induce a small perturbation in the value 
of <p. Assuming that, as r — * oo, cf> — > <pf,, where <fib is the value of (f> in the background, then the outer approximation 
is defined by the assumption — <f>b)/(j)b\ < l/\n + 1|. We may therefore write: 

/ M 3 \ " +1 / Af 3 \ " +1 

-nXM 3 f — J ~ -n\M 3 l—j + m 2 b (cp - </> b ) + O ((</>/& - l) 2 ) 



M X n+2 



where 

ml = Xn(n+ l)M' z I ^ 

V <p 

is the mass of the chameleon in the background. The assumption that |(</> — 4>b)/<f>b\ < + 1| is essentially the 
same assumption as was made in the linear approximation in section UlI Al In the linear approximation, however, this 
assumption was required to hold all the way up to r = 0. All that is required for the pseudo-linear approximation to 
work, is that the outer approximation be valid for all r > r out , where r out is any value of r less than j-j n . This is to say 
that, we need there to be some intermediate region where both the inner and outer approximations are simultaneously 
valid. 



Outside of the body (j> obeys: 



dr 2 r dr \ <p J M p i 

For the outer approximation to remain valid as {r — > oo, (j) — > </>;,} we need 

„ A m^Y' +1 = ^. 



Solving for cj> in the outer approximation, we find <fi ~ <fi* where: 

A p —m\,r 

<t>* = <t> b -— . (A3) 

r 

A is an unknown constant of integration. It will be determined through the matching procedure. 



3. Matching Procedure 



We assume that there exists an intermediate region, r out < r < Ti n , where both the inner and outer approximations 
are valid. This region does not need to be very large. All that is truly needed is for there to exist an open set, about 
some point r — d, where both approximations are valid. We shall consider what is required for such an open set 
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to exist in section lA~4l below. For the moment we shall assume that it does exist, and evaluate (j>o and A. In the 
intermediate region we must have 

cf, ~ ~ cf,* , 

by the uniqueness of asymptotic expansions. The requirement that S(r) <C 1 ensures that m^r <C 1 in any intermediate 
region. Expanding <jf to leading order in vn^r and equating it to <j) we find that 




A = 

3M p i ' 

Now that the previously unknown constants of integration, A and <j>Q, have been found, we can evaluate the conditions 
under which an intermediate region actually exists. 



4. Conditions for Matching 

For the inner approximation to be valid we must certainly require that (cj>(r = 0) / <pc)~^ n+1 ^ <C 1. This is equivalent 



to: 



(m c i?) 2 «2|n+l|((p c / /9b ) 1 /(»+i)-l) 



(A4) 



It is also interesting to note what is required for the pseudo-linear approximation to be valid outside the body i.e. 
without requiring it to be valid for r < R. This gives the weaker condition: 

(m c i?) 2 «3|n+l||(p c /p h ) 1 /(" +1 )-l 

If this weaker condition also fails then, irrespective of what occurs inside the body, we cannot even have a solution 
where (j> ~ cf> outside of the body. As we show below, the condition that a body have a thin-shell is equivalent to the 
requirement that this weaker condition fail to hold. 

For an intermediate matching region to exist, there must, at the very least, exist some d such that, in an open set 
about r = d, both the inner and outer approximations are valid. We must also require that the inner approximation 
be valid for all r < d, and that the outer approximation hold for all r > d. 

For the outer approximation to hold for all r > d we need: 



(n + !)(</>* -0 6 ) 



< 1. 



Using A = (3p c /3M p i and eq. (|A3p we can see that this is equivalent to: 



(m c i?) 2 5 « ( - 



Pb 



l/(n+l) 



(A5) 



Using this condition, we define d m i n to be the value of d for which the left hand side and right hand side of the above 
expression are equal: (m c R) 2 R/3d min — (pb/ Pc)^ 1 ^^ 1 ^ ■ For the outer approximation to be valid at d we require 
d > d 7n i n . 

For the inner approximation, cj> ~ <f>, to hold in the intermediate region we require that, for all R < r < d: 

n+l 



R 3 



> 



dr' 



dr"r" 



(j)(r") 



We must also have that (4> c /4>(r — 0)) n+1 i.e. that condition (|A4|) holds. Provided that this is the case then, for all r 
in (R, d) , we need 



R :i 



> 



dr' 



dr"r" 



n+l 



(A6) 



If both eqs. (|A4[) and (|A6|) hold then the inner approximation will be valid for all r < d. We evaluate eq. (|A6[) 
approximately as follows: For n < —4 we define r = d' by <f> — 4>b = 4>b- For all r in (d',d) we approximate <f> in the 
above integral by and for all r in (R, d') we approximate </> by (</> — fa). 

When n > 0, eq. (|A4[) implies that 4>t > (4> — 4>b) ^ C and so such no d' exists. When n > and condition (|A4|) 
holds, we can therefore find a good estimate for the validity of the pseudo- linear regime by setting tf> = <fo in the above 
integral. We consider the cases n < —4, n = —4 and n > separately below. 
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a. Case: n < —4 



We shall deal with the subcases d! < R and d! > R separately. 
Subcase: d! < R 

The subcase where a" < R includes those circumstances where the linear regime is valid (see section MI Al) . Since 
we already found, in section MI Al how <p behaves when linear approximation holds, we shall only consider what 
occurs when the linear approximation fails. If the linear approximation fails, then the outer approximation must 
break down outside the body i.e. d m i n /R > 1. 

Evaluating eq. (|A6[) . we find that we must, at the very least, require that: 

£<1, 



R 



which, for all n, is equivalent to: 



m c R < V3 ( — 



C R < y/3\n+l\ 



l/2|n+l| 

El 
pb 



The second criteria is just the statement that d! < R. When p b <C p c , this latter condition is more restrictive than 
the former. 

Subcase: d! > R 

From the definition of d! we find d'/R = (m c R) 2 /(3\n + 1|) (/0 c /p 6 ) 1/|n+l1 and d! = d min /\n + It follows that 
d > d m i n > d! . When d! > R, an intermediate region will exist so long as: 

3d 3 p b , 3 ^(m CJ R) 2 V" +1| <<1 



2R 3 p c (n + 4)(n + 3) V3|n+1 

The smaller d is, the more likely it is that this condition will be satisfied; we therefore evaluate the condition at 
d = d m i n . Both of the terms on the left hand side are positive, and so for an intermediate region to exist we must 
have: 

/„ \ (n+4)/6(n+l) / X (n+4)/6(n+l) 



(n + 4)(n + 3)V /2 '" +11 



m c R < y/3\n + 1| 

For n < —4, the second of these conditions is usually the more restrictive. However, since n < —4 implies 
n < —6, this second condition is itself, in general, less restrictive than requiring either (<fi c /cf>(r = 0)) n+1 -C 1 or 
{<t> c /4>{r = R)) n+1 « 1. 

Conditions 

Putting together all of the conditions found above, we see that, for the pseudo-linear approximation to be valid all 
the way from r = oo to r = for n < —4, we must, at the very least, have: 

/ /m \ ("+4)/3(n+2) 1/2 \ 

m c R <mm[(18)W^J , ^TT\ | (m b /m c ) 2 /l"+ 2 l - l| , 

where we have used p c / Pb — {4>b/4>c) nJrl — (m e / 'mb) 2t " n+1 ^' / (™+ 2 ). 
Thin-Shell Condition 

If the above condition fails to hold, then it is a sign that non-linear effects have become important. However, even 
when non-linear effects are important, we only expect them to be so in a region very close to the surface of the body 
itself. In section MI CI we assumed that our body had a thin-shell and considered the behaviour of the field close to 
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the surface of the body. We found there to be pronounced non-linear behaviour near the surface of such bodies. This 
implies that the pseudo-linear approximation must break down for r > R for bodies with thin-shells. We further 
found that the assumption that the body had a thin-shell required: 



n + l 



T./2+1 



m c R > 1. 



For a body to have a thin-shell we must therefore require that this condition hold and that the pseudo-linear approx- 
imation would break down for r > R: these are the thin-shell conditions. It is almost always the case that the latter 
of these two conditions is the more restrictive. In section UlIBI we found this condition to be: 



(m c R) > min (18) 



1/6 



m b 



(n+4)/3(n+2) 



l/|n+2| 



l/2\ 



Since this condition implies m c R 3> 1, it is the thin-shell condition for n < theories. 



b. Case n > 

The case of n > is actually slightly simpler than the n < —4 one because we cannot have d! > R. Other than 
that, the analysis proceeds in much that same way as it does for the n < —4 case. For this reason, we will not repeat 
the details of the calculations here. 



Conditions 

For the pseudo-linear approximation to be valid all the way up to r = we must, at the very least, have: 

^+4 



ni. B mm [ i 18) 1/6 ( ^ 



6(n+l) 



\/2(n + 1) (Pc/pb)^ 1 ~ 1 



1/2' 



When pb <C p c , the most restrictive bound comes from the second term on the right hand side. 
Thin-Shell Condition 

As in n < theories, a body with have a thin-shell provided that m c R 3> 1 and non-linear effects are important 
near the surface of the body which implies that the pseudo-linear approximation breaks down outside the body. This 
latter condition in fact implies the former and is therefore the thin-shell condition for n > theories. This condition 
reads: 



\m b J 



(n+4)/6(n+2) 1/2 \ 

,V3(n + l) KW 1/( " +2L l 



m c R > min 

where the second term on the right hand side is usually the more restrictive when p b <C p c 



c. Case n = —4 



The n = — 4 case, i.e. <j) A theory, requires a more involved analysis. The reason for this is that, unlike the n < — 4 
theories, there does not exist a solution to this theory where cj) ~ — A/r as r — > oo. If we propose such a leading 
order behavior for the inner approximation, then it can be easily checked that the next-to-leading order term dies 
off as ln(r)/r, i.e. more slowly than the leading order one. This means that, for some finite r, the next-to-leading 
order term will dominate over the leading order one. When this happens the inner approximation will break down. 
It can also be checked that higher order terms will always die off more slowly than the terms of lower order. This 
complication will only manifest itself, however, when the conditions for the pseudo-linear approximation fail, or 
almost fail, to hold i.e when (m c R) is large. 

Inner approximation and matching 

To avoid these difficulties, we shall use a different form for the leading order behaviour for <fi when n = —4. We 
write: 

4> (j) b a{r)e- mbr 
4>c 4>c m c r 
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To leading order in the inner approximation we neglect terms of C(mf,r) and smaller. The field equation for <f> is then 
found to be equivalent to: 

d 2 a da a 3 
dy 2 dy 3 

where y = yo + m(r/i?). The outer approximation is still given by <f> ~ <j>* . To perform the matching we need to know 
the large r behaviour of the inner approximation, and the small r behaviour of the outer one. Solving for a we find 
that: 



Hv) ( ^ ' 



The next-to-leading correction to a is: 



2y 

It is clear that, as we would wish, the next-to-leading order term dies off faster than the leading order one. We shall 
see below that, for this approximation to be valid near the body, we need yo 0.634. This approximation also breaks 
down when m^r 3> 1. Matching the inner and outer approximations in the same manner as we did before, we find for 
r>R 

j>_ _ <h + a{mm(r,l/m b ))e- m » r 
4>c 4>c m c r 

Inside the body, r < R, we assume that, at leading order, 4>/<fi c behaves in the same way as if did for all the other 
values of n: 



b (m c R) 2 r 



C C 6(n + l)R 2 
By requiring <f> to be C± continuous at r = R, we find yo to be 

{m c Rf nr i / r~3 



2y 3 W 2y 

and so 

<h _ <h {m c R) 2 1_ 

4>c 4>c 18 m c R V 2y ' 



(A8) 



Conditions 

For this approximation to be a valid approximation, we must firstly require that the next-to-leading order 
correction to a is always small compared to the leading order one. This implies 

yo > 0.508, m c i?< 3.13. 

We must also require <fi(r = 0) = <j>o *C <ftc which gives (assuming Pb/p c *C 1): 

yo> 0.634, m c i?< 2.915. 

This is the stronger of the two bounds. 

When the field equations are solved numerically we find that the form for cf> given above is an accurate approxima- 
tion whenever m c R < 1. 

Critical Far Field 

The form of the far field for bodies with thin-shells in n = —4 theories is examined in appendix [5] below. However, 
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it is interesting to note that, even in the pseudo-linear approximation for n — —4, there is already the first hint of 
/3-independent behaviour in the far field. When m^r -C 1, we found that 

( 2/0 + ln(r/i?))- 1 /2 e ^^ 
9 ~ 9b 



2V2Ar 

Thus, when hi(r/R) yo (provided m^r -c 1), we have, to leading order: 



2 v /2Aln(r/i?)r' 

which is manifestly independent of C and hence also of (3. In n ^ —4 theories, /3-independent critical behaviour was 
reserved for bodies with thin-shells. When n — —4, we can see that leading order, /3-independent behaviour in the far 
field [r 3> R) can occur for all yo i.e. all values of (m c R). However if (m c R) <C 1, this critical behaviour will only be 
seen for exponentially large values of r/R. If, however, m c R > 2 then <f> will be /3-independent at leading order for all 
r/R> 10. 



APPENDIX B: FAR FIELD IN n = -4 THEORY FOR A BODY WITH THIN-SHELL 

Using eq. (|22[) and the other results of section lTlI CI we find that, if a body has a thin-shell, then for < (r — R)/R <C 
1 we have 

1 



-V2A(r-i?) + 4-' 

when n = —4. In appendix [Al we saw that, far from the body and when m^r <C 1, we have: 

, (y + \n{r/R))- 1 ' 2 e- m » r 



2V2Ar 

We can find the value that yo takes for a thin-shelled body, and hence determine the behaviour of the far field, by 
matching the leading order large r behaviour of the first expression to the leading order behaviour of the second 
expression as r — > R. This gives: 

Ayo = 1 => y = 1/4. 

Numerical simulations show that this tends to be a slight over-estimate of the true far-field. We find therefore that 
far from a body with thin-shell, the chameleon field is approximately given by: 

(1 + 41n(r/i?))- 1 / 2 e- m " r 
9~ 9b 



2Ar 

Indeed, since we are far from the body, it is almost always that case that r/R ^> 1.3 and so 41n(r/i?) 3> 1. As a 
result, it is usually a very good approximation to take 



2 N /2Aln(r/i?)r 

We note that, whenever the pseudo-linear approximation holds, yo = 1/4 is equivalent, by eq. (|A8|) . to an effective 
value of m c R of (m c R) e ff — 4.04. Numerical simulations confirm that there is pronounced, thin-shell behaviour 
whenever m c R > 4. When n — —4, the condition for a body to have a thin-shell is therefore m c R > 4. This condition 
is shown to be a sufficient condition for thin-shell behaviour in section MI CI above. 

It is clear that, far from the body, </> is independent of both /3 and the mass of the body, M.. There is a very weak, 
log-type, dependence on the radius of the body, R. The strongest parameter dependence is on A. We can use the 
form of <fi given above to define an effective coupling, /3 e //. Peff is defined so that 

f3 eff (r/R)Me- m >> r 



AnMpir 
It follows that, for r > 2R: 
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APPENDIX C: EFFECTIVE MACROSCOPIC THEORY 

The analysis of the effective field theory for a macroscopic body proceeds along the same lines as the analysis that 
was performed to find the field about an isolated body in section HTT1 and appendices [Al and [Bl A detailed discussion 
of precisely what is meant by 'effective', or 'averaged', macroscopic theory is given section lTVl above. For the purposes 
of this appendix, our aim is to find the average value of the chameleon mass inside a body in a thin-shell. 

We firstly consider what is required for linearisation of the field equations to be a good approximation. We then 
introduce a pseudo-linear approximation. Finally, we consider what we expect to see when non-linear effects are 
strong. We assume that the macroscopic body has a thin-shell, and is composed of spherical particles of mass m p 
and radius R. The average inter-particle separation is taken to be 2D. The average density of the macroscopic body 
is: p c — 3m p / AnD 3 . The density of the particles is p v — 3m p / AnR 3 . We label the average (i.e. volume averaged over 
a scale > D) value of <fi deep inside the body by < (j> >. We shall assume that the effect of the other particles, on a 
particle at r — 0, is sub-leading when r <C D. In general, the surfaces on which d<p/dr — will not be spherical, but 
their shape will depend on how the different particles are packed together, ft will, however, make the calculation much 
simpler, and easier to follow, if we assume that the field about every particle is approximately spherically symmetric 
for < r < D, and that at r = D, d<fi/dr = 0. We take everything to be approximately symmetric about r = D. 
Wc define <fi(r = D) = <p c : m c = m^ c ). We argue below that < <fi >« <fi c , and that the average chameleon mass is 
approximately equal to m c — m<^(<^ c ). Our aim therefore is to find <p c and m c . 

We shall see below that it is possible, for all values of the parameters {D, R}, to construct an outer approximation 
that is valid near r = D. The outer approximation will be valid so long as (<fi — 4> c )/4>c <C l/|n + 1|. In the linear 
regime, this outer approximation will be valid everywhere. In the pseudo-linear and non-linear regimes, however, the 
outer approximation will only be valid for r > d m i n where 

mlD 3 



When n ^ —4, we shall find that R <C D implies that we always have m c D <C I. When n = —4, R <C D implies 
that m c R < is always true. It follows that (d m i n /D) 3 <C 1 • The volume in which the outer approximation holds 
is: 

Vout = MD Z - d 3 min )/3. 

Since (d m i n /D) 3 <C 1, V ou t ~ 4irD 3 /3 i.e. the entire volume of the region < r < D. This means that the volume 
averaged value of <fi, and (</>), will be dominated by the value <f> takes in the outer expansion. Since <fi w <p c in the 
outer expansion, and ~ m c , it follows that < (f> >~ 4> c and < m c . 

Throughout this appendix, we shall therefore refer to (j) c as the average value of </>, and m c as the average chameleon 
mass. The averaged behaviour of cf>, in a body with a thin-shell, is entirely determined by (f> c and m c . Our aim of 
finding an effective macroscopic theory is therefore equivalent to calculating cf> c and m c . 

If linear theory holds inside the body then (f) c — where: 



^(lin) = M f 



nXM p iM 3 



We also define: 

(f> P = M 



p Pp y 1/( " +1) 

nXM p iM 3 J 



1. Linear Regime 

We write 4> = <f> c + (f>\, where \<pi/<p c \ <C I , and <j>(r — D) — cj> c . Linearising the equations about </>o, one obtains: 

+ - + mc(0l/0c) + ^m H{R - r) ' (C1) 

and 



3(3 M b m 2 c D 3 



n+l 



4irR 3 4> c (n + 1)R 3 \ J 1 
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For this linearisation of the potential to be valid we need, just as we did in section IlIIAl |</>i/</> c | <C l/|n + 1|. Solving 
these equations is straightforward and we find that for r > R: 

d>i cosh(m c (r — D))D sinh(m c (r — D)) 

( n + !)— = 1 • 

<pc r m c r 

Inside the particle, r < R, we have: 

/ \ 71+1 
D 3 ( <t>c \ ( smh(m c r)R 



R3 I \ sinh(m c i?)r 

(cosh(m c (i? — D))m c D + smh(m c (R — D))) sinh(m c r) 
smh(m c R)m c r 

For d<fi/dr to be continuous at R = d, we need: 

n+l 

R 6 m c D cosh(m c D) — sinh(m c D) (C2) 
D 3 m c Rcosh(m c R) — sinh(m c i?) 

The largest value of \(f>x/<j) c \ occurs when r = 0. For this linearisation to be valid everywhere we must therefore 
require: 

m c D cosh(m c _D) — sinh(m c _D) f m c R 




m c R cosh(m c i?) — sinh(m c i?) \ sinh(m c i?) 
(cosh(w c (i? - D))m c D + sinh(m c (fl - D))) 
smh(TO c KJ 

We shall see below that m c R -C 1. Given that m c 7? is small, the left hand side of the above condition becomes: 



m 3 c D 3 
2m r R 



3(sinh(m c £>) — m c D cosh(m c Z))) 



i 3 D 3 



1 + smh(m c D)m c D — cosh(m c D) + 0(m c R) 



For this quantity to be small compared with 1, we need both m c D <C 1, which implies m c R <§; 1, and m^D 3 /2R <C 1. 
For all n ^ —4 we define: 



i + 2 



"• = <»<» + < C4 » 
n(n + 1) \ ^ / 3j3m p \ 



^* " {-Jir-J ) ' (C5) 

we note that -D*/D c = (D c /R)( n+1 ^ 3 . With these definitions the requirement that (m c D) 2 < 1 is equivalent to 
D ^s> D c . When n < —4 we need fl C D* for m^D 3 /R <C 1, whereas when n > we need £) 3> -D* for the same 
condition to hold. Therefore, when n > we need £) ^> £> c , D*, whereas for n < —4 we require D c <C -D <C D*. 

When n > 0, no matter what value R takes, there will always be some range of D for which the linear approximation 
holds. If n < —4, however, we must require that -D* 3> D c for there to exist any value of D for which the linear 
approximation is valid. It follows that, for the linear approximation to be valid for any D (when n < —4) we need 
R > D c , i.e. m^{(j) p )R < 1. 

When n = —4 we need both D <C -D* and: 



fl2 W2 A i/2 ( 3 ff TO P 



< 1. 



This second condition implies both that m c D <IC 1 and Tn^^tpp 

We conclude that, for large enough particle separations, it is always possible to find some region where the linear 
approximation holds when n > 0. However, when n < —4 we must also require that m < f,(<p p )R <C 1 for there to be 
any value of D for which the linear approximation is appropriate. 

Whenever the linear approximation holds, it follows from m c D -C 1 and eq. (|C2[) that we must have: 
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2. Pseudo-Linear Regime 

The pseudo-linear approximation proceeds in much the same way as it did in section fill Bl and appendix [XJ for an 
isolated body. Near each of the particles we can use the same inner approximation as was used for a single body. 
This is because, near any one particle, the other particles are sufficiently far away that their effect is very much 
sub-leading order. 

Inner and Outer Approximations 

The inner approximation (for n ^ —4) is therefore: 

= A - 

where A is to be determined by the matching procedure. We will deal with the n — —4 case separately later. 

Previously, the outer approximation was defined so that it remained valid as r — > oo. In this model, however, we 
are assuming that everything is symmetric about r = D, and so we need only to require that the outer approximation 
remain valid up to r = D. Near r — D, we assume that <fi ss <p c and linearise about <p c . Requiring that d<j)/dr = 
when r = D, we find: 

,| 1 / cosh(m c (Y — D))D sinh(m c (r — D)) 




n + 1 \ r m c r 

The outer approximation, as defined above, is also good for n —4. 

Matching 

For the pseudo-linear approximation to work we need, just as we did in appendix ^ there to exist an intermediate 
region where both the inner and outer approximations are simultaneously valid. We will discuss what conditions this 
requirement imposes shortly, but, before we do, we shall assume that such a region does exist, and match the inner 
and outer approximations. Matching the inner and outer approximations we find: 

1 H — (1 - cosh( m c D) + m c D sinh(m c L>)) , (C6) 

1 + n 

3(cosh(m c D)m c D — sinh(m c D)) 




Conditions for matching 

For the outer approximation to remain valid in the intermediate region, where r ss d say, we require that, for all 
r G (d,D): 



1 — cosh(m c Z?) + m c D sinh(m c _D) 



ID 3 3(cosh(m c d)m c D — sinh(TO c D)) 



M 



{m c Df 



< 1. 



This is equivalent to m c D <C 1 



. We must also require 



ID 3 



3d 



< 1. 



We define d m i n to be the smallest value of d for which the above condition holds. 

For the inner approximation, ~ (j), to hold in the intermediate region, we require that conditions, similar to those 
that were found in the isolated body case, hold (see section Ull Bl and appendix Specifically, we require that for 
all r in (R, d) 



Y > / <//•' / <//■"/■" 



ra+l 



(C8) 



We also require that (0 c /0(r = 0)) ,l+1 <C 1. We note that ((j)c / 4> P ) n+1 — P P /p c — D 3 /R 3 . We consider the subcases 
n < —4, n > and n = —4 separately. 
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a. Case n < —4 



The analysis proceeds in the same way as it did for an isolated body in appendix [X] We find that we must require 

/ R \3/\n+l\ s 



mlR 2 = { ^ 



1+4 
i + l 



< 2\n+ II 



If we only wish for the pseudo-linear approximation to remain valid up to r = R, then we must require the weaker 
condition </>(r = R)/(j> < 1. This is equivalent to 



j n 2 / D c 
"' = ~~R 



< 3k- 



"-(5) 



3/|n+l 



Provided that either of these conditions hold, the conditions for the outer approximation to be valid in the intermediate 
region are automatically satisfied. 

Whenever the first of the above requirements holds, the pseudo-linear approximation will give accurate results. 

When the latter (and weaker) of the two conditions fails, we expect pronounced non-linear behaviour near the 
particles. When this happens the far field induced by each particle becomes /?-independent. We will discuss how this 
affects the values of (f> c and m c in section [C3l below. 



m c D <C 1 implies 



via eq. IC71 It follows that 



therefore, looks precisely like it did in the linear regime. 



(0( Zm )). The resulting macroscopic theory, 



b. Case n > 

The analysis for the n > case proceeds in the much same way as it did for a single body (see appendix For 
the outer approximation to hold, we require that: 



D/D c > 1, 



which implies m c D <C 1. We also need 



~~D 



< 



2(77+1) 1 



R\—^ 
D. 



If we only wish to require that the pseudo-linear approximation hold up to r = i?, then we can relax this second 
condition to: 



~~D 



< 



3(77+1) 1- 



Whcn the first of the two conditions holds, the pseudo-linear approximation gives accurate results. Whenever the 
second condition fails, we expect pronounced non-linear behaviour near the particles. When this happens, we expect 
that the far field (in r>R) will attain a critical form. We discuss the consequences of this in section [C"3l below. 



As in the n < —4 case, when the pseudo-linear approximation holds we have m c D <C 1 and so 



This 



implies m c 



lin) 



c. Case n — —4 

The case of n = —4 is, as it was in the one particle case, the most complicated to study. However, the analysis 
proceeds almost entirely along the same lines as it did in the one particle case. We find, using the results of appendix 
\K[ that the field near the particles will behave like: 



A 



75- R 

m c r 



(C9) 
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where y — yo + ln(r/i2), and, provided non-linear effects are small i.e. yo > 0.6, we find yo to be given by: 

9 V 2t/ 3 ^22/0, 

When n = —4, we have that m^(<p lm )D = m < f > (4> p )R. 

Using this, and matching the inner approximation to the outer one, we find: 

A = H — (1 - cosh(m c D) + m c Dsmh(m c D)) , (CIO) 

1 + n 



2(y +\n(D/R)) 



= (cosh(m c D)m c D - sinh(m c £>)) /3. (Cll) 



where y(D) = y + \n(D/R). 

Whenever non-linear effects arc small we have yo > 0.6, which implies 

iY/3< nr 

2yo) V 2yo 

and so 

"3~ (m4^ p )R) 2 _ {m^ 1 ^ D)f 



2y 9 

Since y > 0.6 implies J ^ < f .6 and so it follows eq. (|C11|) that: 



(m c Df » (m (^ im ^)) E 



f + \n{D/R)/y ■ 



The above approximation is very accurate when it predicts m c D < 1 , and even when m c -D > 1 it gives a good estimate 
for m c D. 

For macroscopic, everyday, bodies with densities of the order of 1 — 10gcm~ 3 , we tend to find \n(D/R) w 11. If 
linear theory is to give a good estimate of m c , we need: 

ln{D/R)/y «l^-t, »ll=> m c (<j> (lm) ) < 1.5. 

More generally, linear theory gives a good approximation to m c whenever: 

2(m 4> ( ( f> ilm ^Df ln(X>/iZ)/243 « 1. 

Therefore, unless D/R is improbably large (D > Re 243 / 2 ), we will have we will have m c ps mi 1 " 1 ^ whenever mi hn ^ D < 
1. When m c (</>( /m )) > 1.5, we actually move into the a regime of /3-independent, critical behaviour, more about which 
shall be said below. 



d. Summary 

When the pseudo-linear approximation holds (and m^{(j^ lln ^)D < 1), we have found that m c « m^((/)' m ). As a 
result, despite the fact that there does not exist an everywhere valid linearisation of the field equations, linear theory 
actually gives the correct value of m c , at least to a good approximation. 



3. Non-linear Regime 



When the pseudo-linear approximation fails it is because non-linear effects have become important near the surface 
of the particles, and they have developed thin-shells of their own. Far from the particles we expect that the field will 
take its critical form as given by eqs. (f2"4"k -c). We can find the value of (f> c , in this case, by matching our the outer 
approximation for <p to the critical form of the far field around the particles. 
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When n < —4 we will have critical behaviour whenever: 



„ + 4 

R 



> 3|n+ II 1 



3/|n+l| N 



When this happens the far-field (r 3> i?, r <C I?) behaviour of the chameleon takes its critical form. We have found 
this to be well approximated by: 



cf)~ (f) = Acj) c 



7(n) 



l/\n+2\ 



(MR) 
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Performing the matching to the outer approximation we find that: 



A = 1 — — (1 — cosh(m c D) + m c D s'mh(m c D)) , 



7(n) 



l/|n+2 



We therefore have that: 



(MR)^ 



(n + l)m ( 



(cosh(m c D)m c Z) — sinh(m c .D)) 



n(MDf ( M\ n+1 3 (cosh(m c D)m c D - smh(m c D)) 



(m c D)3 



n(MD) 3 (M 



n+l 



cent) „ v3F+ti m ™ f 7W\ ™ 



D 



When n > 0, a similar analysis finds that 



D 



> 



3(n + l) 1- 



R — 
D 



3 



and that 



We therefore find 



i(MDf ( M\ 



n+l 



MR 



n(n + 1)M 2 



2\ V("+2) 



(crit) 



y/3|n+l| 



1/2 



In the n = —4 case, critical behaviour will actually emerge whenever ln(D/R)/yo > 1. Non-linear effects are still 
responsible for this critical behaviour but it is not necessarily the case that the particle have developed thin-shells. 
Indeed, the thin-shell condition requires that m < / > (<f> p )R > 4, whereas critical behaviour (for ln(D/i?) w 11 as is 
typical), emerges whenever m$(4> p )R > 1.4. When the particles do have thin-shells, yo = 1/4 (see appendix |B|) and 
so critical behaviour is seen whenever D > 1.3i?. 

We find the (3- independent critical mass for n — —4 theories using eq. (jCllj) and taking yo <C \n(D/R). It follows 
that, the critical value of m c for n = —4 theories is given by: 



m 



(crit) 



X/D, 



where X satisfies 



X cosh X — sinh X 



3^3 



When In(D/iJ) = 11, we find m^ rit D w 1.4. 
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4. Summary 

In this appendix, we have performed a very detailed analysis of the way in which the chameleon behaves inside a 
body that has a thin-shell and which is made up of many small particles. In this way, we have been able to calculate 
the average chameleon mass inside such a body. Despite the in depth nature of the analysis, our results can be 
summarised in a very succinct fashion. 

The average mass, m c , of the chameleon field inside a thin-shelled body of average density p c is 

/ — „l U>. \ T 

m r — mm 



where R is the radius of the particles that make up the body, and 2D is the average inter-particle separation. The 
critical mass, m[ mt ', is given by: 



in , 



Jo-it) ra x j D ^ n = -4 i 
where q(n) = mm((n + 4)/(n + 1), 1), S(n) = 1 if n > 0, S(n) = ( 7 (n)/3) 1 / 2 l" +1 l if n < 0, and 

3V3 



X cosh X — sinh X 



The dependence of rric D vs. ln(D/i?) is shown in FIG. [4] We can clearly see that D ^> R implies mi cnt ^ D <C 1 
when n ^ —4. When n — —4, m'f™^ can be seen to be less than V3 whenever ln(D/i?) > 2. 

APPENDIX D: EVALUATION OF q FOR WHITE-DWARFS 

In this appendix, we evaluate for white-dwarfs under a more accurate approximation than that used in section 
IVIII This accurate evaluation allows for the effect of a chameleon on the general relativistic stability of white dwarfs 
to be studied in more detail. The effect of the chameleon on compact objects such as these is considered in section 

Em 

The chameleon contribution to the energy of the white-dwarf was found, in section IVIII to be: 

w 71+1 f JS ^ 

W$ = / d 3 x—p. 

n J Mpi 

For the chameleon theory to be valid, we must require that the chameleon force is weak compared to gravity inside the 
white dwarf. We shall therefore ignore the chameleon corrections in the behaviour of p when evaluating W<p, because 
they must be sub-leading order. In this way, we are able to accurately find the leading order behaviour of W^. By 
leading order we mean: ignoring general relativistic and chameleon corrections. It is important to stress that we are 
not ignoring special relativistic corrections, which are very important; we are merely assuming that the only volume 
force inside the white dwarf is, to leading order, Newtonian gravity. This approximation is similar to the one used in 
[13] to evaluate the general relativistic corrections to the energy of the white dwarf. 

We shall assume that the polytropic equation of state, P = Kp T , holds everywhere. We will also approximate the 
white dwarf to be spherically symmetric. Defining p — p c 9 p and r — a£, hydrostatic equilibrium implies that, to 
leading order, 6 satisfies the Lane-Emden equation [57] : 

i_d_ 2 de 



where 



(p+l)KM^p, 

47T 



2 Jl/P-I) 



1/2 
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TABLE II: Values of for different n 



n 


«0 




u 




-12 


0.8493 




4 


1.523 


-10 


0.8207 




5 


1.410 


-8 


0.7786 




6 


1.337 


-6 


0.7110 




7 


1.286 


-4 


0.5853 




8 


1.249 


1 


3.613 




9 


1.220 


2 


2.144 




10 


1.197 


3 


1.720 




11 


1.178 



We define p c to be the density of matter in the centre of the star. The boundary conditions 9(0) — 1, 6'(0) = 1 follow. 
The index p is related to T by T = 1 + l/p. This equation must be solved numerically. The surface of the star is at 
r = R = a£i, which is defined to be the point where p = 0. We will mostly be interested in the case of relativistic 
matter, r = 4/3 — > p = 3. For p — 3 we have £i = 6.89685. The chameleon potential is then given by: 

n + 1 f3<f>(p c ) 3 



n M p i J Q 

We evaluate this integral numerically assuming a relativistic equation of state (i.e p = 3) and find: 

71+1 MPc)., 
= a 0^T7 Abator 

n Mpi 

where M star = m u N is the mass of the star. The values of a$ are given in table |TTJ As n — » ±oo we find — » 1. 
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